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Radio  wave  propagation  in  the  earth's  crust  is  reviewed.  This  has  been 
studied  by  many  investigators  who  usually  assume  the  crust  to  be  a planar 
slab  of  low  loss-tangent,  bounded  on  top  and  bottom  by  regions  of  large  loss- 
tangent.  Available  conceptual  depth  profiles  of  complex  dielectric  constant 
of  the  lithospheric  waveguide  are  reviewed  and  theories  extant  summarized. 
Sources  were  usually  VED's  inserted  through  vertical  bore-holes.  For  ocean 
emplaced  antennas  located  near  the  ocean-lithosphere  boundary,  the  sources 
are  horizontally  polarized  and  propagation  is  predominantly  by  means  of  lateral 
waves.  Calculations  are  reported  for  a wide  range  of  conditions.  The  useful 
characteristics  of  lateral  wave  transmission  are  discussed.  The  measured 
properties  of  scaled  waveguide  propagation  that  have  been  reported  could  not 
accurately  scale  the  losses  in  the  overburden  since  it  was  represented  by  a 
metal  plate.  The  principles  of  electrodynamic  similitude  are  delineated  and 
two  different  experiments  using  similitude  with  different  depth  profiles  are 
discussed,  t 

To  study  the  effects  of  boundaries  on  transmission  (including  wall  losses) 
and  simultaneously  to  make  measurements  of  embedding  and  moving  antennas  in 
them,  the  laboratory  models  should  have  "walls"  composed  of  a liquid  with 
suitable  excitation.  For  an  orderly  build-up  of  modelling  measurements , a 
laboratory  set-up  is  described.  Included  would  be  measurements  with  boundaries 
that  are  sloping,  undulating  or  irregular.  Though  available  materials  may 
limit  accurate  modelling  according  to  similitude  scaling,  meaningful  and 
illuminating  information  could  be  obtained  for  better  understanding  of 
lithospheric  transmission.  The  more  complex  depth  profiles  require  a much 
more  elaborate  set-up,  patterned  possibly  after  that  of  lizuka. 
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1.  Introduction 

Interest  in  the  electrical  and  geophysical  properties  of  earth's  crust 
has  been  recently  encouraged  by  the  possible  existence  of  a lithospheric 
plane  waveguide  with  conductivities  lower  than  previously  accepted.  The 
region  of  the  earth's  crust  forming  this  low  conductivity  waveguide  was 
assumed  to  be  a planar  slab,  10  to  30  km  thick,  bounded  on  top  by  an  over- 
burden (or  ocean)  of  high  conductivity  and  on  the  bottom  ("underburden")  by 
a high  conductivity  mantle  or  Moho  (as  it  is  often  referred  to  in  the 
literature) . A review  of  the  crustal  properties  was  given  at  an  ONR 
Lithospheric  Workshop  in  mid-March,  1978,  chaired  by  J.  Heacock  with  papers 
presented  by  G.  Keller,  G.  Simmons,  P.  R.  Bannister,  and  C.  Cox  amongst 
others. 

There  have  been  some  excellent  reviews  of  the  field  in  recent  years. 

For  example.  Wait  (1971)  in  the  AGU  Monograph  No.  14,  covered  principally 
the  waveguide  analytical  approaches.  And  3annister  (1978)  gave  an  excellent 
review  of  the  idealized  waveguide  (homogeneous  profile)  of  Wait  and  Spies 
plus  the  results  of  Field  and  Dore  for  non-homogeneous  profiles  of  Levin, 
Keller,  and  Housley.  We  wish  merely  to  supplement  here  some  of  their 
writings. 

While  most  of  the  analytical  work  tended  to  treat  waveguide  modal 
solutions,  there  exists  the  possibility  of  lateral  wave  transmission  between 
ocean  emplaced  HED  transmitting  and  receiving  antennas  located  near  the 
ocean-lithosphere  boundary. 

It  was  desired  to  study  the  feasibility,  in  laboratory  models,  of 


measuring  the  properties  of  lithospheric  transmission,  including  wall  losses, 
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using  the  principles  of  electrodynamic  similitude  for  scaling.  Scaled 
models  have  been  used  for  a homogeneous  waveguide  in  one  case,  and  in  another 
for  a non-homogeneous  profile  of  conductivity  increasing  with  depth.  Both 
cases  used  metallic  upper  plate  upper  walls  which  do  not  accurately  scale 
the  properties  of  the  overburden  or  ocean.  There  is  a question  of  available 
materials  for  such  modelling. 

In  Chapter  2,  we  present  separately  the  various  available  conceptual 
depth  profiles,  starting  from  the  1961  work  of  Wheeler  and  that  of  Mott  and 
Biggs  in  1963,  both  works  being  expositions  of  concepts  of  communication 
from  land  via  "A  deep  waveguide"  in  the  crust  to  deeply  submerged  submarines. 
There  follows  a series  of  descriptions  of  several  non-homogenecus  depth 
profiles. 

Chapter  3 contains  a summary  review  of  the  various  theoretical  works  on 
the  profiles  of  Chapter  2.  Most  of  these  concern  modal  solutions  for  wave- 
guides with  plane  boundaries;  results  are  often  illustrated  with  boundaries 
having  infinite  conductivities.  In  some  cases,  the  lower  boundary  is  sharp; 
in  others  the  complex  dielectric  constant  is  tapered  to  increase  with  in- 
creasing depth.  Finally,  the  work  with  antennas  emplaced  in  the  ocean  near 
the  ocean-lithosphere  boundary  is  cited  by  way  of  introducing  the  next 
section. 

In  Chapter  4 is  an  expose  with  numerical  calculations  of  lateral  waves 
from  an  HED  source  emplaced  in  the  water  near  the  ocean-lithosphere  boundary. 

A previous  development  of  King  and  Sandler  for  an  HED  in  water  below  air, 
which  extended  the  work  of  Banos,  is  extended  to  the  case  where  the  lithosphere 
replaces  the  air.  A wide  range  of  parameters  were  used  for  calculating  the 


fields 


Chapter  5 contains  a brief  development  of  expressions  for  electrodynamic 
similitude.  Examples  are  given  for  laboratory  frequency  and  distance  scaling 
by  Brown  and  Gangi  for  a homogeneous  depth  profile  of  dielectric  constant  and 
conductivity,  plus  some  results  of  lizuka  with  the  non-homogeneous  profiles 
achieved  with  an  agar-agar  gel  diffused  with  sodium  chloride. 

In  Chapter  6 is  contained  a discussion  of  lateral  wave  modelling  in 
combination  with  measurement  of  losses  in  the  walls  for  this  and  for  wave- 
guide excitation.  A suggested  laboratory  set-up  could  be  used  for  sloping 
boundaries  or  boundaries  that  are  rough  or  undulating. 

A sunmary  and  conclusion  with  recommendations  are  contained  in 
Chapter  7. 

There  follows  the  list  of  references  and  supplemental  bibliography  in 
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2.  Lithospheric  Profiles 


2.1  Introduction 

Most  of  the  earlier  lithospheric  waveguide  depth  profiles  were  based 
on  the  laboratory  measurements  of  samples  of  rock  as  functions  of  tempera- 
ture and  pressure.  With  few  exceptions,  emphasis  was  placed  on  the  depth 
variation  of  conductivity  o;  this  was  due  probably  to  the  belief  that  the 
expected  exponential  attenuation  would  be  low  enough  only  at  low  frequencies 
where  the  loss  tangent  p = 60  c XQ/er  is  quite  large  (a  in  Si/m,  \Q  the  free 
space  wavelength  in  meters,  and  er  the  relative  dielectric  constant) . In 
analyses  of  these  profiles,  most  investigators  had  assumed  er  to  be  constant 
with  depth  in  the  lithospheric  wave guide,  with  values  of  4 to  10. 

Earlier  analyses  assumed  the  waveguide  to  have  a depth  of  10  to  40  km. 
On  top,  it  was  bounded  by  an  overburden  or  sea  of  such  high  loss-tangent 
that  this  region  was  represented  by  a highly  conducting  metallic  plate  in 
many  analyses  and  models. 

At  the  bottom,  the  earlier  analyses  assumed  the  conductivity  towards 
the  hot  "Moho"  discontinuity  to  be  so  high  that  it  again  was  represented  by 
a highly  conducting  plate.  These  earlier  analyses  assumed  a homogeneous 
parallel  plate  slab,  somewhat  idealized,  to  calculate  the  fields  in  the 
guide.  Refinements  grew  as  more  complicated  profiles  of  the  crust  developed. 
We  trace  here  the  profiles  from  the  early  one  of  Wheeler  to  that  of  Wait, 
then  to  Levin  and  coworkers,  to  Keller,  to  Housley  and  to  the  work  of 


Simmons. 
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2.2  Wheeler’s  Profile 

Wheeler's  (1961)  early  paper,  interestingly  enough,  concerned  itself 
with  a means  of  communicating  to  a submarine  near  an  ocean  floor  from  a 
shore  based  source,  A sketch  of  the  concept  is  shown  in  Figure  2.1  (a); 
the  transmitter,  shown  in  Figure  2.1  (b) , passes  through  about  1 km  of 
overburden  and  then  extends  2 km  in  a borehole  into  the  waveguide  where 
it  excites  a TEM  mode.  The  conductivity  profile  is  shown  in  Figure  2.1  (c) 
as  a function  of  depth  and  temperature.  The  boundaries  are  roughly  2 and 
20  km  deep,  where  the  rock  conductivity  may  be  10“6  to  lO”^1  Si/m,  a useful 
value  being  10-8  Si/m.  The  relative  dielectric  constant  was  6.  The  "upper 
wall",  perhaps  1 to  several  kms  thick  has  a range  of  4 Si/m  for  sea  water 
down  to  10-4  Si/m  in  rather  dry  minerals.  The  lower  region  was  termed  the 
"inverted  ionosphere"  or  "thermal  ionosphere"  by  Wheeler  (1961)  with  tem- 
peratures of  300°  to  600°  C. 

2.3  Mott  and  Biggs'  Profile 

A couple  of  years  later,  Mott  and  Biggs  (1963)  also  considered  VLF 
propagation  below  the  sea.  Their  concept  is  shown  in  Figure  2.2  (a)  for 
the  coiranunicating  region,  the  upper  boundary  1 km  thick  being  the  wet, 
highly  conducting  continental  surface  or  the  ocean.  The  propagating  region 
was  assumed  to  consist  of  dry,  low-temperature  low  conductivity  granites  or 
basaltic  rocks,  20  to  70  km  thick.  The  bottom  was  formed  by  ultra-basic 
rooks,  having  higher  conductivity  because  of  the  high  temperature. 

Propagation  takes  place  between  2 vertical  electric  dipoles  (VED) 
extending  1 or  2 km  into  the  low  conductivity  region,  with  an  upper  value 
of  10“6  Si/m  for  long-distance  communication. 
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(a)  Conmunicaticn  frcm  buried  \TD 
source  S via  deeo  vaveouide 
under  the  ocean  to  receiver  R 
located  in  the  sea. 


(fc)  FVetch  of  TTD  source  antenna  S 
inserted  into  the  vave«-uide. 
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(c)  Profile  of  oonductivitv  and 
tenoerature  with  derth. 
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Fio.  2.1  Vheeler's  concent  and  conductivitv  orc^ile  for  ccrmunicati on  under  the  ocean 
via  a deeo  vaveouide  to  a submerged  receiver.  (From  Vbeeler,  1961) . 
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Earth  conductivity  estimates  are  shown  in  Figure  2.2  (b) . The  curve 
labeled  number  "1"  indicates  measured  earth  conductivities  of  10-3  for 
poorly  conducting  earth  to  4 Si/m  for  sea  water.  Curve  "2"  is  from  Hughes 
(1953)  for  olivines,  with  McDonald's  (1959)  estimates  of  temperature 
distribution.  Curve  "3"  is  an  estimate  of  maximum  conductivity,  by  Garland 
and  Webster  (1960)  for  western  Canada  pre-Cambrian  basement  rocks.  Curve  "4" 
is  that  of  Lahiri  and  Price  (1939)  while  curve  number  "5"  is  due  to  Cantwell 
and  Madden  (1960) . Some  of  these  curves  are  included  in  a composite  one 
drawn  by  Iizuka  (1968) , discussed  later  on.  The  final  conductivity-depth 
profile  assumed  by  Mott  and  Biggs  (1963)  is  shown  in  Figure  2.2  (c) . In 
essence  this  is  a two-region  model,  the  bottom  being  10  - 20  km  deep  with 
a bottom  region  having  a conductivity  10“^  to  10-7  Si/m,  and  a dielectric 
constant  of  6 with  a top  region  of  sea  water  with  4 Si/m  for  conductivity, 
perhaps  1 to  2 km  thick. 

2.4  The  Profile  of  Brown  and  Gangi 

Brown  and  Gangi  (1962,  1963)  and  Gangi  (1966)  gave  a diagram  of  earth 
conductivity  as  a function  of  frequency  and  depth,  as  shewn  in  Figure  2.3. 

The  surface  and  near  surface  overburden  had  conductivities  ranging 
from  10"1  to  10”4  Si/m,  with  an  average  of  10”2  to  10-3  Si/m.  The  relative 
dielectric  constant  had  a wide  range  of  values,  with  a typical  value  of  4 to 
10  being  widely  used. 

For  the  crust,  the  values  were  inferred  from  tables  in  Von  Hippel 
(1954)  for  the  effects  of  temperature  and  water  content  of  quartz  and  sili- 
cates (pages  303-314  of  his  book) . The  low  conductivity  value  of  granite 
is  frequency  sensitive  and  the  authors  show  values  for  1,  10  and  100  kHz. 


At  30  km  depth,  the  conductivity  increases  rapidly,  temperature  being  about 
300°  C,  relatively  independent  of  frequency.  Near  the  Moho  the  conductivity 
is  about  4 x 10-3  Si/m.  Relative  dielectric  constant  was  of  the  order  of  5, 
down  to  30  km. 

In  their  theoretical  studies,  the  "real  earth"  situation  was  approximated 
by  a uniform  slab  between  two  highly  conducting  planes  (See  next  Chapter) . 

The  conductivities  chosen  were  according  to  frequency  in  the  low  conductivity 
region. 

2.5  Levin's  Profile 

An  idealized  cross-section  of  the  earth's  crust  was  delineated  and  shown 
by  Levin  (1966)  in  Figure  2.4  illustrating  the  crust  under  a continent  vs.  an 
ocean  basin,  as  well  as  the  continental  shelf.  Under  the  continent  is  granite 
underlain  by  basalt  down  to  the  Moho.  Proceeding  through  the  shelf,  the 
granitic  rock  disappears  leaving  basalt  under  the  ocean  and  volcanoes.  The 
Moho  may  be  35  km  deep  under  the  continent  and  about  10  km  under  the  ocean. 

After  a considerable  discussion..  Levin  (1966)  showed  a trial  profile  of 
both  conductivity  and  relative  dielectric  constant  with  depth.  Levin  (1966) 
cited  the  measurements  of  rock  characteristics  by  Acker  and  Mueller  (1966) 
from  3 km  drill  holes  associated  with  a full  scale  propagation  experiment. 

His  conductivity  profile  is  summarized  in  Levin  (1971)  and  shown  in 
Figure  2.5.  The  variation  of  dielectric  constant  is  discussed  later  in 
the  paper  by  Schwering,  Peterson,  and  Levin  (1968) , and  is  overlooked  by 
many  authors. 

In  his  AGARD  paper,  Levin  (1966)  concluded  that  "full  scale  studies  of 
lithospheric  propagation  must  be  carried  out  in  boreholes  nearly  10  km  deep 
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to  assure  access  to  the  dielectric  zone  (loss-tangent  less  than  unity)  of 
the  crust  over  the  entire  length  of  the  path". 

In  Figure  2.5,  two  interpretations  are  shown  for  conductivity  vs. 
depth.  The  curve  labeled  "dry"  is  more  nearly  in  favor  in  view  of  recent 
consideration  of  laboratory  measurements  and  interpretations,  especially 
those  regarding  the  influence  of  microcracks  (Simmons,  1978). 

2.6  The  Schwering,  Peterson,  and  Levin  Model  - Inverse  Square  Profile 

Schwering,  Peterson,  and  Levin  (1968)  analyzed  lithospheric  propaga- 
tion based  largely  on  the  Levin  (1966)  model,  the  dry  portion  of  conductivity 
of  Levin  (1971)  shown  in  Figure  2.5.  The  variation  of  dielectric  constant 
and  conductivity  are  shown  in  Figure  2.6  as  solid  lines  as  the  "experimental 
values" ; the  dotted  curves  are  used  for  their  analysis,  described  in  the  next 
Chapter. 

It  was  convenient  to  distinguish  three  zones  in  the  crust,  with  elec- 
trical constants  given  for  1 kHz,  and  shown  in  Table  2.1,  as  well  as  other 
postulated  characteristics.  The  waveguide  is  Zone  2 in  essence. 

2.7  Wait's  Exponential  Profile 

After  studying  the  performance  of  electrical  constants  in  a number  of 
ways.  Wait  (1966c) , offered  a unique  profile  suitable  for  analysis  (Wait 
1966a,  1971),  to  attempt  a quantitative  result.  His  model  is  illustrated 
in  Figure  2.7.  In  a cylindrical  coordinate  system,  the  plane  z = 0 sepa- 
rates the  highly  conducting  overburden  having  a complex  relative  dielectric 
constant  Kg,  the  overburden  top  being  at  z = -h.  The  basement  has  a rela- 
tive dielectric  constant  K (z)  which  is  complex,  having  a constant  relative 


Figure  2.6  Sch"erircr,  Petersen  arc  Tevir  096P)  profiles  cf  relate  dielectric 
constant  and  conduct  ivi  tv.  Fee  also  Fait,l?71.  ’"'ashed  cT.,jr”e  is  tue 
rodel  orofile;  solid  curve  is  that  used  Ccr  theoretical  analvsi*. 
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seJ  n rntarv  .gneous  and  metamorphu.  ' • . pen  wsicpi  grounj- 

v».jic'  Naturjied.  climate  and  terrain  permitting  ('onductivitv  and 
dielectric  permittmiv  controlled  b\  :nc  water  distribution  hence 
characterized  h\  great  .ariabiiitv  high  !o  >m  \aiuev  large  negative 
gradient  witn  depth,  a ■*  ]()' 1 to  10  * .nm  m.  r = JO''  to  10*  . loss 
tangent  ! Probable  depth  extent  about  km 

Zone  2 ! *. f .ilJie > Ham-  jnd  Jr\  Zone  of  tight  g'.innu  basement  oc^ 
more  homogeneous . rc  atoei)  closed  -vstem.  little  or  no  mo'mi^r 
H-O.  : v pica i of  silicate  dielectric  Li  * to  r. -derate  <nlO  ' u>  10  “ 
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dicnts  with  depth  . S<  I Zone  exiend<  to  14*  or  4<»  km  depth 

Zone*  L*>*tr\  />i  but  not  Bottom  /one  of  the  w .ot.  p,*s\iblv  extendme 
nto  the  mantle  Anhydrous  «iiicate  m:^crj  .*w.nh  age  vf  p utonK 
tav.es  Intrinsic  jnd  ionic  charge  .omerv  rad  to  moderate  to  nigh 
c mducUMtv  and  permittivity  «j>I0  * r?  I(H)|.  with  »en  Urge 
positive  gradients  «$  :>  1 

V.»/<  The  \ aiues of  c . and  ^ applv  to  a trequenev  ot  1 >Hz 


Table  2.1  Descriptions  of  crustal  zones  according  to  Schwering,  Petersen, 
and  Levin  given  in  their  1968  article,  Schwering,  Petersen,  and 
Levin  (1968) 
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dielectric  constant  and  a conductivity  which  increases  exponentially  with 
depth.  The  reference  depth  was  chosen  for  convenience  as  that  depth  zd 
where  the  loss  tangent  was  two;  for  VZF  it  appears  that  z^  is  of  the  order 
of  30  km.  At  very  great  depths,  where  the  loss  tangent  greatly  exceeds 
unity,  Y.(z)  varies  exponentially  with  depth.  Near  z * 0,  the  ideal  case 
would  have  K(0)  real. 

The  results  of  using  this  profile  will  be  given  in  the  next  Chapter 
where  he  compares  his  results  with  those  of  other  authors, 

2.8  Iizuka's  Profiles 

Iizuka  performed  a number  of  experiments  using  a "typical"  lithospheric 
profile  deduced  from  a number  of  sources.  The  typical  case  of  a profile  of 
loss-tangent  with  constant  relative  dielectric  constant  is  shown  in  Chapter 
5,  Figure  5.6.  for  a dielectric  constant  of  4 and  a typical  frequency  of 
5 kHz.  The  compilation  of  conductivity  profiles  of  Iizuka  (1968  ) is  shown 
in  Figure  2.8,  They  were  deduced  from  Wait  (1963) , Wheeler  (1961) , Mott 
and  Biggs  (1963),  Keller  (1963),  Ames,  Frazier,  and  Orange  (1963),  Cantwell 
and  Madden  (1960) , Watt,  Mathews,  and  Maxwell  (1963) , deBettencourt  and 
Frazier  (1963),  Garland  and  Webster  (1960)  among  others. 

In  Figure  2.8,  the  conductivity  is  plotted  as  a function  of  depth  to 
depths  of  40  km  or  slightly  more.  The  conductivity  ranges  from  10-6  to  10-1 
Si/m,  the  larger  values  being  near  the  earth's  surface  or  at  great  depths. 

Since  he  preferred  an  analysis  based  upon  loss-tangent,  he  noted  "lower 
limit  of  tan  5"  near  a conductivity  of  10“8  Si/m  and  a higher  limit  near 
10  km  depth  with  a conductivity  of  about  10“4  Si/m.  (For  his  "typical" 
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Figure  2.P  TizuFa's  simarv  of  conductivity  profiles  fmrr  various  sources , ci rca  ]Qf3.  TizuFa  (196P  ).  See  text. 
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loss  tangent  curve,  Iizuka  (1968  5 assumed  a conductivity  of  slightly  less 
than  10  7 Si/m  with  a dielectric  constant  of  4,  at  a frequency  of  5 kHz) 

2.9  Keller's  Profile 

Keller  (1963,  1971,  1978a, b)  has  studied  in  depth  the  electrical  crustal 
properties.  His  most  recent  profile  of  conductivity  is  shown  in  Figure  2.9. 

It  has  been  widely  used  and  quoted.  It  appeared  first  in  Gallawa  and  Haidle 
(1972),  as  a private  communication,  again  in  Field  and  Dore  (1973),  discussed 
by  Keller  (1978a)  at  the  ONR  Lithospheric  Workshop  and  is  reprinted  here  from 
Keller  (1978b) . 

While  the  value  of  10“4  Si/m  at  very  shallow  depths  appears  low  (Field 
and  Dore,  1973;  Bannister,  1978),  the  curve  down  to  about  10  km  is  based  on 
Appalachian  well-log  data  and  an  extrapolation  thereof  towards  minimum  values 
of  the  order  of  10”7  Si/m  at  10  km.  The  lower  portion  has  increasing  values 
from  those  of  granite  to  underlying  basalt.  A rather  complete  discussion  is 
given  in  Field  and  Dore  (1973) . 

If  water  saturated  microcracks  exist  in  rocks  at  5 to  8 km  depth,  mini- 
mum conductivities  must  be  greater  than  10-5  to  10-6  Si/m,  and  practical  long 
range  communication  does  not  appear  feasible.  However,  recent  data  concerning 
microcracks  by  Simmons  (1978)  tend  to  indicate  conductivities  could  be  lower, 
say  10  8 to  10  Si/m,  than  that  previously  believed  (10_/  Si/m). 

2.10  Housley's  Profiles 

Housley  (1973)  made  laboratory  measurements  on  the  conductivity  of  dry 
rock,  see  Simmons  (1978).  Results  show  even  lower  conductivities  than  those 
of  Figure  2.5  of  Levin  and  Figure  2.9  of  Keller.  Using  the  new  conductivity 
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Figure  2.9  Keller’s  profile  of  conductivity  and  temperature  vs.  depth. 

From  Field  and  Dore  (1973) , first  as  a private  communication 
in  Gallawa  and  Haidle  (1972),  and  discussed  by  Keller  (1978a, b) 
See  also  Bannister  (1978) . 
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data,  the  profiles  are  similar  to  those  of  Keller  in  Figure  2.9,  for  depth 

less  than  10  km,  and  are  shown  in  Figure  2.10.  At  greater  depths,  much 

_8 

lower  conductivities  are  indicated  (less  than  10  Si/m) . 

The  three  profiles  correspond  to  geotherms  in  three  different  heat 
flow  provinces  (i.e.,  BR  = Basin  and  Range;  EUS  = Eastern  United  States; 

SN  = Sierra  Nevada) . The  thermal  gradients  for  the  Keller  model  was  16°/km, 
while  for  the  Housley  provinces  they  were  BR  (24°/km) , EUS  (14°/km) , and 
SN  (9°/km) , where  the  temperature  is  in  (°C) . 

2 .11  Comment 

There  are  many  differing  profiles  of  conductivity  vs.  depth  which  are 
not  necessarily  definitive.  Continued  experimental  work  remains  to  be  done 
in  the  laboratory  - see  Simmons  (1978)  - and  in  the  field  - see  Keller  (1978a) . 


3.  Summary  Review  of  Theories 


3.1  Introduction 

When  viewing  the  number  of  more  recent  depth  profiles  of  electrical 
constants  in  Chapter  2,  such  as  those  of  Keller  in  Figure  2.9  and  of  Housley 
in  Figure  2.10,  interest  has  been  revived  in  the  possibility  of  long-distance 
communication  via  a sub-surface  waveguide.  The  minimum  conductivities  are 
several  orders  of  magnitude  lower  than  those  postulated  in  the  early  1960 
era.  The  lower  'onductivity  values  have  received  support  from  studies  of 
the  effects  of  microcracks  by  Simmons  (1978). 

A number  of  review  papers  (Rjanzcev  and  Shabel 'nikov,  1966;  deBetten- 
court,  1966;  Levin,  1966,  1971;  Wait,  1966c,  1971;  Gabillard,  1966;  Gabillard, 
Degauque,  and  Wait,  1971;  and  Bannister,  1978)  have  been  published  on  the 
various  mechanisms  dealing  with  lithospheric  propagation,  including  ground 
waves,  "up-over-and-down, " and  sub-surface  waveguide  mechanisms  plus  com- 
munication parameters  of  antennas,  noise,  etc.  Here  we  shall  briefly  sum- 
marize those  papers  dealing  principally  with  the  crustal  waveguide  mechanism. 
Many  focus  attention  on  attenuation  of  the  modes,  a few  give  examples  of 
field  strength  variations. 

3.2  Some  Early  Works 

It  is  interesting  to  realize  that  two  of  the  earliest  papers  considered 
the  possibility  of  using  the  lower-loss  region  of  the  crust  for  the  purpose 
of  illustrating  a potential  mechanism  for  communication  with  submarines  at 


great  depth 


Wheeler's  (1961)  very  early  concept  was  shown  in  Figure  2.1.  He 

imagined  the  existence  of  a very  deep  waveguide.  With  a relative  dielec- 
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trie  constant  of  6 and  a conductivity  of  10  to  10  Si/m,  he  obtained 
a total  attenuation  of  the  TEM  launched  wave  of  0.067  dB/km. 


Mott  and  Biggs  (1963)  also  considered  undersea  communication  between 

two  vertical  electrical  dipoles  (VED)  both  inserted  into  lower  conductivity 

crustal  rock.  They  envisioned  propagation  at  frequencies  of  the  order  of 
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100  kHz  to  take  place  in  a plane  half  space  of  conductivity  10  to  10  Si/m 
bounded  on  top  by  the  layer  (ocean)  having  a conductivity  of  4 Si/m.  They  go 
on  to  assume  the  lower  region  attenuates  the  wave  but  does  not  affect  the 
propagation  mechanism  - surface  wave  along  the  interface  like  normal  surface 
wave  propagation  in  air  above  the  sea  but  where  the  air  has  been  replaced  by 
the  rock.  See  also  Carolan  and  deBettencourt  (1963) . At  lower  frequencies, 
less  than  10  kHz  for  example,  the  lower  region  behaves  as  a waveguide,  and 
"that  a TEM  waveguide  mode  will  be  important",  as  Wheeler  (1961)  had 
postulated.  Their  conductivity  profile  is  shown  in  Figure  2.2(c). 


3.3  Slabs  with  Uniform  Depth  Profiles 


For  analytical  studies  of  their  depth  profile  shewn  m Figure  2.3, 
Brown  and  Gangi  (1962,  1963)  and  Gangi  (1966)  assumed  the  crustal  waveguide 
to  be  homogeneous  slab  bounded  by  two  half-spaces  which  were  highly 
conducting.*  For  their  analysis,  the  boundaries  were  assumed  to  have 
effectively  infinite  conductivity.  The  waveguide  was  excited  by  a vertical 
electric  dipole. A mode  series  representation  was  used  to  obtain  the  field. 


‘Reference  may  be  made  to  Figure  5.1  for  a sketch  of  the  guide 


r 
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Calculations  were  made  of  the  vertical  electric  field  strength  versus 
range  for  several  low  frequencies  for  a guide  depth  H = 35  km,  a relative 
dielectric  constant  of  4 and  zero  antenna  depths.  A typical  result  is 
shown  in  Figure  3.1  for  1,  3,  and  10  kHz,  the  conductivity  values  corres- 
ponding to  their  profile  shown  in  Figure  2.3.  For  this  configuration,  the 
modal  cut-off  frequency 


f (m)  = m v2/2H  (3.1) 

where  v^  is  the  phase  velocity  in  the  crustal  waveguide  and  where  m is  the 
mode  number,  whence  f (m)  = 2.1m,  kHz.  Hence,  for  1 kHz,  only  the  TEM 
mode  (m  = 0)  propagates  readily;  for  3 kHz  there  are  two  modes  and  for  10 
kHz,  there  are  five  modes  (m  = 0,  1,  2,  3,  and  4),  propagating  with  low 
attenuation.  Interferences  occur  because  each  mode  has  a different  attenua- 
tion and  phase  shift. 

Using  Wait's  (1957)  mode  theory  for  the  earth-ionosphere  waveguide 
with  lossy  ground,  but  replacing  the  air  by  a low-loss  waveguide,  deBetten- 
court  (1962,  1966)  calculated  the  vertical  electric  field  versus  range  for 
infinitely  conducting  walls  for  a guide  width  of  20  km,  relative  dielectric 
constant  of  4,  and  for  conductivities  of  0 and  10  6 Si/m;  the  dipole  moment 
was  16,900  amp.m.  (A  plot  of  the  results  may  be  seen  in  Figure  6 of  Wait, 
1966c)  . 

Spies  and  Wait  (1971)  and  Wait  and  Spies  (1971)  have  extensively  con- 
sidered the  wave  in  an  idealized  waveguide.  From  a general  solution  of 
multiple  horizontal  stratified  layers  (Wait,  1962) , Spies  and  Wait  (1971) 
chose  a restricted  3 region  case,  shown  in  Figure  3.2,  in  which  the 
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"overburden"  region  had  infinite  conductivity,  the  crustal  "waveguide" 
(region  2)  had  uniform  conductivity  of  c2  = a and  relative  dielectric 
constant  K2,  and  the  "sub-basement"  (region  3)  had  a conductivity  O3  = ac 
and  relative  dielectric  constant  of  Kc,  and  was  infinitively  thick.  The 
waveguide  has  a thickness  h2  = h = H.  The  boundaries  of  the  various 
regions  are  sharp.  In  their  calculations  using  a Newton  method  of  suc- 
cessive approximations,  Spies  and  Wait  (1971)  and  Wait  and  Spies  (1971) 
considered  frequencies  from  1 to  100  kHz;  these  were  expended  later  for 
frequencies  down  to  1 Hz  (Wait  and  Spies,  1972). 

Wait  (1966,  1971)  demonstrated  the  principle  of  "super-position  of 
losses,"  which  also  was  used  by  Watt,  Leydorf,  and  Smith  (1966).  For 
modes  of  sufficiently  low  attenuation,  the  total  exponential  attenuation 
of  a mode  can  be  considered  as  the  sum  of  three  terms,  i.e»,  for  a mode  m 

At  (m)  = Ap  + + Aj  (3.2) 

where  the  three  components  are  respectively  the  plane  wave,  upper  boundary, 
and  lower  boundary  contributions  to  the  total  attenuation  rate. 

We  are  generally  concerned  with  a crustal  waveguide  having  a low- 
loss-tangent  region  whence  it  can  be  shown  that 

Ap  = 1.64  x 106  n2/JT2  dB/Kn  (3.3) 

and  is  independent  of  frequency  if  the  electrical  "constants"  are  constant. 
For  a given  plane  slab  boundary,  the  wall  attenuation  contribution  depends 
upon  the  guide  loss-tangent  and  whether  the  boundary  is  electrically  thick 
or  thin.  Expressions  based  on  approximation  similar  to  those  of  Wait  (1963) 
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have  been  given  for  several  cases  by  Watt  et  al.  (1966)  . As  aforementioned, 
we  are  generally  interested  in  the  crustal  waveguide  which  has  a low  loss- 
tangent,  and  the  boundary  layers  have  a high  loss-tangent,  and  are  electri- 
cally thick.  For  this  case,  if  o is  the  conductivity  of  one  of  the  uniform 
slab  boundary  layers,  then  Watt  et  al  (1966)  give  the  boundary  attenuation 
for  the  TEM  (m  = 0)  mode  as 

A0  (0)  = 2.29  x 10“2  (K2f/a  )Vh  dB/Km  (3.4) 

Whilst  for  the  TMj  Mode  (m  = 1) 

-‘h 

Ab  (1)  = 4.58  x 10"2  (X2f/°Q)h  [l-(  fc/f)2J  /H  dB/Km  (3'5> 

where  f is  f (m  = 1)  from  equation  (3.1) 
c c 

In  Figures  3.3  and  3.4,  taken  from  Spies  and  Wait  (1971),  the  total 
attenuation  for  modes  of  lowest  attenuation  are  shown  as  functions  of 
frequency  for  a guide  thickness  of  H = 20  km,  and  upper  layer  of  infinite 
conductivity,  the  waveguide  with  a dielectric  constant  K = 9 and  conduc- 
tivities of  10-6  or  10~2  Si/m,  and  the  sub-basement  layer  with  a uniform 

dielectric  constant  of  K =15  and  a conductivity  of  10“2  Si/m  (Fiaure  3.3) 

c 

or  10“3  Si/m  (Figure  3.4). 

In  a manner  similar  to  that  of  Bannister  (1978) , we  made  calculations 
for  a frequency  of  10  kHz,  the  model  cut-off  frequency  being  2.5m  kHz. 

(The  loss-tangent  of  the  waveguide  region  is  very  small.  The  plane  wave 
attenuation  A^  for  = 10“^  Si/m,  = 4 works  out  to  be  .0546  dB/Km, 
whilst  for  0 = 10"6,  it  becomes  0.546  dB/Km.  The  plane  wave  attenuation 

is  critically  dependent  on  conductivity. 


Frequency  (kHz) 


Fig. 3.4  Sane  as  Fig. 3. 3,  but  for  (T.  » 10-3  Si/m  (^pies  and  '’ait,l°71) 
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For  the  Spies  and  Wait  (1971)  model  calculations  for  Figure  3.2,  the 

upper  layer  conductivity  is  infinite  so  A (m)  = 0 and  A (m)  = A ♦ As(m). 

u c pi 

For  H =*  20  km,  Kj  = 9,  for  the  mode  (m  = 1)  , f » 2.5  kHz.  For 
= 10 ~2  Si/m  and  f = 10  kHz,  the  boundary  attenuation  is  .007  dB/Kta 
while  for  ag  ~ ^2  = Si/m  it  is  .023  dB/Km.  The  higher  the  conductivity , 

the  lower  the  attenuation  or  leakage  through  the  wall. 

Similar  calculations  were  made  for  the  TEM  (m  = 0)  mode.  Some  of  the 
results  are  shown  tabulated  in  Table  3.1  where  the  values  in  parentheses 
are  rounded  off  tabular  values  taken  from  Spies  and  Wait  (1971) ; the  approxi- 
mate attenuation  AT  (0)  and  A^,  (1)  can  be  compared  with  Figures  3.3  and  3.4; 
the  agreement  is  quite  good  here. 

It  is  interesting  to  consider  one  aspect  of  the  attenuation  of  the 
lowest  order  mode,  the  TEX  mode  (m  = 0) , with  reference  to  equations  (3.2), 
(3.3),  and  (3.4).  If  the  homogeneous  slab  of  relative  dielectric  constant 
K.,  and  conductivity  C.,  has  a thickness  of  2h_  = 2h  * 2H,  but  is  bounded  on 
top  and  bottom  by  homogeneous  half-spaces  of  the  same  conductivity  oa  and 
dielectric  constant  the  total  attenuation  (0)  will  be  unchanged 

from  that  derived  for  the  configuration  of  Figure  3.1.  This  will  be  referred 
to  in  Chapter  6 in  the  discussion  of  model  measurements  of  wall  losses  with 
an  image  plane  for  region  1.  This  point  has  also  been  noted  by  Wait  and 
Spies  (1972)  in  their  discussion  of  ELF  crustal  waveguide  attenuations  for 
the  lowest  order  mode. 

The  model  just  discussed  is  mindful  of  the  analysis  for  a three  layer 
structure,  shown  in  Figure  3.5(a),  as  considered  by  Viggh  (1963).  His  slab 
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TAPLF  3.1 


Aoproximate  versus  exact  modal  attenuations  for  modes  of  lew  attenuation,  for 

for  idealized  wavecruide  shewn  in  Fiaure  3.2 


n-0*- 

V>  h2 

■ H = 20  km.  f = 10  kHz 

€r2  = *2  = 9 

er3 

= K3  = 15 

TEM  (ir*=0) 

TV]_ 

(ir=l) 

c r 

3 = 

10"2  Fi/m 

2 

A3(0)  At(0) 

^2 

PP 

*3(1) 

Aij>(1) 

10-6 

.546 

.004  .55  (.55) 

10-6 

.546 

.007 

.55  (.57) 

10“7 

.055 

.004  .059  (.059) 

10*7 

.055 

.007 

.062  (.064) 

°~3  = 10-3  ci/p 


10-6  .546  .011  .56  ( .56) 

10-7  .055  .011  .066  (.071) 


10~f  .5^6  .023  .57  (.59) 

10-7  .055  .023  .079  (.090) 


For  coroarison,  the  ficures  in  oarentheses  are  from  Fnies  and  T-Tait,l°71,  tables. 
The  attenuations  A are  in  dP/km. 
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of  depth  D with  conductivity  c2  = 10“6  Si/m  was  bounded  by  identical  top  and 
bottom  half  spaces  with  conductivities  = 10“2  Si/m;  all  three  regions 

had  the  same  relative  dielectric  constant  of  10.  The  slab  was  electrically 
thin  so  that  only  the  TEM  mode  was  considered  dominant.  An  example  of  the 
attenuation  as  a function  of  thickness  is  shown  in  Figure  3.5(b)  where  fre- 
quency is  the  parameter;  both  scales  are  logarithmic. 


Waveguides  having  low  loss-tangents  but  bounded  by  layers  on  top  and 
bottom  have  been  considered  by  Gabillard  (1966)  and  in  a review  paper  given 
by  Gabillard,  Degauque,  and  Wait  (1972).  Essentially  these  papers  refer  to 
the  work  in  France  where  they  were  able  to  demonstrate  the  propagation 
effects  in  mine  galleries  and  in  a uniform  gypsum  layer. 


3.4  Non-homogeneous  Depth  Profiles 

3.4.1  Wait's  Exponential  Profile 


Besides  the  aforementioned  works  of  Wheeler  (1961)  and  Mott  and  Biggs 
(1963)  another  early  paper  was  that  of  Wait  (1963,  1971) who  considered  the 
analysis  of  the  exponential  depth  profile  shown  in  Figure  2.7  in  order  to 
obtain  quantitative  results.  For  z positive,  the  complex  relative  dielec- 
tric constant  was  chosen  for  convenience  to  be  exponential  in  form  and 
given  by 


K (z)  = Kq  £l-j2  exp  B(z-z(j)j 


(3.6) 


where  z^  is  the  waveguide  "depth"  which  is  discussed  below. 
For  z = 0,  at  the  top  of  the  waveguide 


Kq  [l-j2  exp  (-Bzd)]  * Kc 


I 


K(0) 


(3.7) 
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provided  Bzd  exceeds  3.  With  zd  of  the  order  of  30  km  and  B of  the  order 
of  0.2/km,  this  restriction  is  true,  i.e.,  the  loss-tangent  is  extremely 
small  and  K0  is  essentially  a real  quantity  at  z = 0. 

For  z = zd, 

K(z)  = K0(l-j2)  = Kq/5  exp  (3.8) 

That  is,  the  loss  tangent  at  z = zd  is  2. 

For  very  large  depths 

K(z)  = -j j K (z)j  = -j2KQ  exp  B(z-zd)  (3.9) 

Since  Kq  is  real,  the  conductivity  may  be  expressed  as  a function  of 
depth  as 

a(z)  = Kq  exp  (z-z^J/30  XQ  (Si/m)  (3.10) 

where  X0  is  the  wavelength  in  free-space  in  meters.  These  expressions 
assume  that  this  real  relative  dielectric  constant  and  conductivity  are 
not  functions  of  frequency. 

At  z = 0,  where  the  loss-tangent  of  the  guide  is  very  small,  we  desig- 
nate the  wavelength  as 

\'  = X /K1*  (3.11)* 

2 o o 

Wait  (1963)  derives  the  vertical  field  strength  for  a VED  in  the  wave- 
guide from  a vertically  directed  Hertz  vector.  From  the  residues  at  the 
poles  of  the  denominator  in  the  expression  for  the  latter,  he  obtains  the 

i 

♦Wait  uses  the  quantity  X0  to  designate  what  we  have  called  X2 
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model  resonance  condition,  finally  expressing  the  field  E as  E W where  E 

z o o 

is  a reference  field.  This  reference  field  is  that  of  a vertical  dipole 
in  a dielectric  half-space  erected  on  a perfectly  conducting  ground  plane 
of  infinite  extent.  The  expression  for  w is  dimensionless  and  contains  all 
the  characteristics  of  the  waveguide  modes.  Because  the  lower  boundary  is 
not  sharply  defined,  the  mode  of  lowest  attenuation  is  not  the  TEM  (m  = 0) 
mode  but  one  of  the  modes  (m  = 1,  2 . . .). 

For  purposes  of  illustration.  Wait  (1963,  1971)  calculated  the  attenua- 
tions for  modes  m = 1 and  2 for  the  special  case  where  the  overburden  was 

so  largely  conducting  that  was  assumed  to  be  infinite.  He  chose  the 

/ 

case  where  - 15  km.  The  attenuations  A^  and  A^  for  the  modes  of  lowest 
attenuation  are  shown  plotted  as  functions  of  the  quantity  B in  Figure  3.6(a) 
and  (b) , respectively.  The  parameter  on  the  curves  is  the  waveguide  depth 

i 

in  "wavelengths"  denoted  by  z^/^. 

Wait  (1963,  1971)  shows  how  these  curves  can  be  scaled  for  other  wave- 
lengths, for  a given  guide  thickness  in  wavelengths.  In  these  articles  and 

I 

another  for  plane  slabs  (Wait,  1966) , he  illustrates  how  the  total  attenua- 
tion for  the  lowest  order  modes  may  be  shown  to  be  the  sum  of  the  three  terms, 
the  plane  wave  attenuation  plus  the  losses  in  the  two  "walls"  as  given  in 
equation  (3.2),  (Wait,  1971).  i 

For  the  dominant  mode  attenuation,  he  notes  the  comparison  of  values  of 
Schwering,  Peterson,  and  Levin  (1968)  for  their  inverse  square  profile  of 
0.256  dB/km  compared  with  that  value  (0.275  dB/km)  they  obtained  using  Wait's 
(1963)  method.  Wait's  (1963)  values,  he  notes,  lead  to  lower  attenuations 
than  those  of  Schwering  et  al  (1968)  because  of  the  steeper  gradients. 

i 
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3.4.2  The  Inverse  Square  Profile  of  Schwering,  Peterson,  and  Levin 

The  depth  profile  of  Schwering,  Peterson,  and  Levin  (1968)  is  shown  in 
Figure  2.6,  the  dashed  curves  being  used  for  analysis.  Both  the  real  rela- 
tive dielectric  constant  K and  conductivity  a vary  as  the  inverse  square  of 
the  depth  z measured  positively  downward.  The  vertical  electric  dipole 
(VED)  is  located  in  the  waveguide  at  z = 0 at  the  upper  wall.  The  upper 
zone  (overburden)  is  assumed  to  be  infinitely  conducting.  The  depth 
dependence  is  such  that  both  K and  a become  infinite  again  at  a depth 
y0  = 33km,  and  the  conductivity  of  the  lower  zone  below  33  km  is  also 
infinite. 

The  complex  relative  dielectric  constant  may  be  written 

K* (z)  = K0 (l-jp)F(y0,z)  . (3.12) 


where  F(ye,,z)  is  the  depth  dependence  factor 


(3.13) 


for  the  depth  z (in  km)  lying  between  0 and  33  km.  The  real  relative 
dielectric  constant  at  zero  depth  is  K = 10  and  the  loss-tangent  is  p = 0.2, 
for  the  values  assumed  by  Schwering  et  al  (1968)  at  a frequency  of  1 kHz. 

For  these  assumptions,  then,  the  depth  dependencies  of  the  real  relative 
dielectric  constant  and  of  the  conductivity  are 


K(z)  = K0  F(y0,  z)  (3.14) 

c(z)  = 1.11  x 10"7  F(y0,  z)  Si/m  (3.15) 

respectively.  If  these  electrical  constants  are  independent  of  frequency, 
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then  the  loss  tangent  at  any  frequency  is 
p(f)  = 0.2/f 

where  f is  the  frequency  in  kHz. 

Schwering  et  al  (1968)  employ  a z-directed  vector  potential  to  obtain 
a wave  equation  for  this  potential  by  using  cylindrical  coordinates  and 
related  expressions  for  Ez,  and  Ef.  The  wave  equation  differs  from  the 
usual  wave  equation  by  the  addition  of  a term  due  to  the  ncr.-homogenecus 
depth  behavior  of  K' (z) . Going  through  the  process  of  separation  of  variables 
for  the  vector  potential,  the  authors  arrive  at  a modal  equation  for  each 
field  component  which  includes  "depth  gain"  functions.  At  the  depth  z = 0, 
the  radial  electric  field  Er  = 0,  and  the  corresponding  expressions  for  E2 
and  H$  are  simplified. 

For  first  three  modes  m = 1,  2,  and  3,  the  attenuation  coefficients  at 
1 kHz  are  calculated  to  be  0.256,  1.17,  and  1.99  dB/km  (taken  from  Table  II 
of  Schwering  et  al,  (1968) , which  lists  values  for  the  first  10  modes) . 

The  magnitude  of  the  vertical  electrical  field  strength  is  shown  plotted 
in  Figure  3.7  for  1 kHz  and  for  the  range  20  to  100  km,  using  the  first  eight 
modes.  The  first  mode  is  dominant  at  distances  in  excess  of  50  km;  the 
effect  of  higher  order  modes  is  hardly  evident  in  the  range  of  20  to  50  km. 

In  Figure  3.8  are  shown  the  frequency  variation  of  the  attenuation  rate 
and  phase  velocity  Vi  of  the  dominant  mode  (m  = 1)  , the  latter  normalized 

O 

to  the  velocity  in  free  space  (3  x 10  m/sec) . There  is  resultant  dispersion 
of  the  signal  waves,  but  it  becomes  small  as  the  frequency  increases. 
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Schwering  et  al  (1968)  considered  the  uniform  depth  profile  solution 

_7 

of  Brown  and  Gangi  (1963) , assuming  K = 10  and  a = 10  Si/m  for  a frequency 
of  1 kHz;  they  obtained  the  attenuation  for  the  least  attenuated  mode  as 
0.052  dB/km.  In  the  SPL  model  u.iis  would  correspond  to  the  dominant  (m  = 1) 
mode  attenuation  of  0.256  dB/km.  This  difference  emphasizes  the  sensitivity 
of  the  propagation  constants  to  the  characteristics  of  the  assumed  depth 
profile.  Of  importance  is  the  effect  of  temperature  gradient;  it  is  zero 
for  the  Brown  and  Gangi  model  and  about  20°C/km  of  depth  for  the  SPL  model. 

3.4.3  Theoretical  Results  of  Field  and  Pore  for  More  Recent  Profiles 

Concerned  with  the  possibility  of  lithospheric  waveguide  communication 
systems  as  a whole.  Field  and  Dore  (1973)  considered  the  fields  propagated 
via  the  more  recent  non-homogeneous  profiles  discussed  in  Chapter  2 plus 
the  effects  of  depth  attenuation  of  atmospheric  noise  and  other  system 
parameters  such  as  antenna  performance,  receiver  noise  figure,  and  data 
rates.  Frequencies  considered  extended  from  100  Hz  to  100  kHz.  Using 
conceptual  profiles  having  dry  basement  rock  indicated  more  favorable  com- 
munication transmission  ranges  in  the  30  - 100  kHz  VLF  band.  Presumably 
this  effect  was  due  largely  to  the  neavier  downward  attenuation  of  atmos- 
pheric noise  at  the  higher  frequencies.  An  excellent  review  of  their  system 
analysis  work  was  given  by  Bannister  (1978). 

We  shall  discuss  only  their  results  for  attenuation  rates,  and  calcu- 
lations of  field  strength  versus  range  and  depth. 

Field  and  Dore  (1973)  used  full-wave  calculations  for  the  inhomogeneous 
conductivity  depth  prefiles.  Pectangular  coordinates  were  employed  for 
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plane  TM  waveguide  modes  propagating  in  the  x-direction  for  a vertical 
electric  dipole  lying  along  the  z-axis.  Expressions  were  developed  for 
the  Ex  and  Ey  electric  field  strength  components  and  a "voltage"  normalized 
magnetic  field  component  H’y  = 120ir'Ky  where  Hy  is  the  usual  magnetic  field 
strength.  While  their  formulations  accounted  for  the  depth  variation  of 
the  real  relative  dielectric  constant,  their  calculations  assumed  a 
constant  value  of  10  for  all  depths  and  profiles. 

Presented  in  Figure  3.9  are  the  calculated  attenuation  rates  for  the 
modes  of  least  attenuation  as  a function  of  frequency  for  the  depth  con- 
ductivity profiles  of  Levin  (1966)  shown  in  Figure  2.5,  of  Keller  (1972) 
in  Figure  2.9,  and  Housley  (1973)  BR  profile  shown  in  Figure  2.10.  The 
least  attenuated  modes  are  TEM  modes  for  frequencies  less  than  about  2 kHz 
and  the  lowest  order  TM  mode  for  higher  frequencies. 

In  Figure  3.9  it  is  seen  that  the  attenuation  depends  markedly  on 
both  the  minimum  conductivity  and  the  gradient  cf  the  depth  profile.  The 
lowest  attenuation  occurs  for  the  Housley  BR  profile,  having  a minimum  con- 
ductivity of  about  5 x 10'9  Si/m.  The  Levin  (wet)  model  with  a minimum  of 
10-5  Si/m  has  the  highest  attenuation.  While  the  Levin  (dry)  and  Keller 
models  have  similar  minimum  conductivities  of  10-7  Si/m,  the  attenuation 
rates  at  1 kHz  are  0.28  dB/km  for  the  Levin  (dry)  model  but  only  0.18  dB/km 
for  the  Keller  profile.  For  the  latter,  the  lower  value  is  due  in  part  to 
the  sharp  transition  in  conductivity  near  a 14  km  depth,  according  to  Field 
and  Dore  (1973) . 

Note  that  for  the  uniform  ideal  waveguide  treated  by  Spies  and  Wait 
(1971)  and  Wait  and  Spies  (1971) , with  a conductivity  of  10“7  Si/m  for  the 
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waveguide,  the  resulting  attenuations  are  2 or  3 times  lower  than  the 
values  shown  in  Figure  3.9  for  the  Levin  (dry)  and  Keller  models.  These 
models  have  more  gradual  "boundaries"  than  the  sharp  boundaries  in  the 
depth  profiles  of  Spies  and  Wait.  If  we  averaged  the  Keller  profile  over 
some  range  of  depth,  an  average  might  be  5 or  6 x 10~7  Si/m;  such  a value 
in  the  Spies  and  Wait  case  would  result  in  a higher  attenuation  value. 

In  Figure  3.10  are  presented  the  normalized  magnitudes  of  the  three 
field  components  versus  depth  for  the  Keller  profile  at  1 kHz.  The  nor- 
malization is  such  that  the  value  of  the  "voltage  normalized"  magnetic 
field  H'y  is  unity  at  the  depth  of  minimum  conductivity  at  10  km.  For  the 
VED  receiver,  the  magnitude  of  E_^  has  a broad  maximum  near  a depth  of  10  km. 

Note  that  Ex  has  a maximum  magnitude  at  a 6 km  depth,  14  db  below  the 
maximum  value  of  E-;  Ex  is  less  than  Ez  at  depths  below  a 4 km  depth  but 
Ex  exceeds  Ez  at  smaller  depths. 

Considering  the  vertical  field  itself,  the  relative  vertical  strength 
versus  frequency  for  several  depths  is  shown  in  Figure  3.11  for  the  Keller 
profile.  The  values  are  normalized  to  that  value  at  10  km  depth  irrespec- 
tive of  the  frequency.  From  these  curves,  one  notes  the  relative  decrease 
in  vertical  electric  field  strength  as  the  receiving  antenna  is  moved 
toward  the  surface,  fcr  a giver,  frequency.  Thus,  at  a frequency  of  1 kHz, 
the  relative  field  strengths  are  13,  32,  and  44  dB  below  that  at  10  km 
depth,  for  depths  of  4,  2,  and  1 km,  respectively.  At  20  kHz,  the  degra- 
dations are  16,  36,  and  52  dB  for  those  same  depths,  respectively.  Field 
and  Core  have  made  similar  calculations  fcr  the  Housley  BR  model.  Fcr 
depths  below  10  km,  the  degradations  in  vertical  electric  field  strength 
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are  nearly  the  same  because  the  two  conductivity  depth  profiles  are  the 
same.  (Of  course,  from  a system  point  of  view,  one  must  take  into  account 
the  depth  variation  of  atmospheric  noise  in  order  to  obtain  the  resulting 
signal-to-noise  ratio  at  the  receiver.  The  results  are  not  very  straight- 
forward and  one  is  referred  to  the  original  report  (Field  and  Dore,  1973) 
for  derivation  and  results. 

3 . 5 Other  Studies  - Ocean  Emplaced  Antennas 

3.5.1  The  Work  of  Raytheon  Co.  (Grossi  et  al) 

Grossi  et  al  at  Raytheon  Co.  (1973)  had  considered  the  fields  from 

transmitting  antennas  emplaced  in  the  ocean  or  in  the  lithosphere  below, 

with  the  receiving  dipole  in  the  ocean  near  the  ocean-lithosphere  interface. 

The  model  employed  is  sketched  in  Figure  6.2,  except  that  the  region  below 

_2 

the  uniform  lithosphere  was  assumed  to  have  a conductivity  of  10  Si/m 
instead  of  the  infinite  value  assumed  by  Frieman  and  Kroll  (1973) . 
Frequencies  of  interest  were  in  the  sub-Hertz  ranee.  An  error  in  the 
original  expressions  for  the  fields  was  corrected.  There  was  very  good 
agreement  for  the  magnetic  fields  between  that  for  a horizontal,  long, 
insulated  antenna  on  the  ocean  floor  (which  antenna  was  termed  a "phantom 
loop"  by  Raytheon)  and  that  for  an  insulated  HED  antenna  when  using  the 
formulations  of  Frieman  and  Kroll  (1973).  The  comparison  was  made  at  a 
frequency  of  0.1  Hz  and  ranges  of  100  and  1000  km;  the  lithospheric  region 
had  a conductivity  cf  10  ^ Si/m. 

3.5.2  The  Work  of  Frieman  and  Kroll 

Frieman  and  Kroll,  1973,  determine  the  modes  and  the  mode  excitation 
in  a flat,  stratified  region  for  both  vertical  and  horizontal  current 


distributions.  The  region  of  interest  lies  between  perfectly  conducting 


plane  boundaries  at  z = -d^  and  z = h where  d^  is  the  depth  of  the  litho- 
sphere and  h is  the  height  of  the  ionosphere,  both  measured  from  the  ocean 

floor.*  The  depth  of  the  ocean  is  dw.  The  parameters  used  are  jw=  4 Si/m, 

-6 

o ^ = 10  Si/m.  With  frequencies  in  the  range  between  f = 0.1  and  10  Hz, 
the  propagation  constants  are,  for  the  lithosphere  Y?  = (jwuc^)15,  for  the 

L 2 1* 

sea  water,  Yw  = (jwuaw)  and  for  the  air,  Ya  = (-<d  ue  ) . By  confining 

their  attention  to  modes  that  propagate  to  distances  large  compared  to  both 
dj  and  h,  Frieman  and  Kroll  require  only  two  modes.  These  are  the  litho- 
spheric mode  which  is  large  in  the  lithosphere  and  the  air  mode  which  is 
large  in  the  air.  Solutions  for  the  magnetic  field  of  these  modes  are  ob- 
tained with  both  vertical  and  horizontal  electric  dipoles  located  on  the 
floor  of  the  ocean.  Their  results  indicate  "that  at  large  distances  the 

HED  is  more  effective  than  the  VED  bv  a factor  (c  /a.)*5  = 2x10^",  They  go 

w ** 

on  to  show  that  the  contribution  from  the  terms  that  involve  the  depth  d^, 
of  the  ocean  are  essentially  the  same  as  for  an  ocean  of  infinite  depth. 
Finally,  in  the  range  where  the  thickness  d^  of  the  lithosphere  is  effec- 
tively infinite,  the  magnetic  field  reduces  to  that  obtained  by  Banos,  1966, 
or  Kraichman,  1970,  for  the  problem  of  two  half-spaces.  This  is  given  by 
the  near-field  formula 
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In  a specific  example,  the  magnetic  field  is 
insulated  antenna  10  km  long  and  4 cm  in  diameter 


calculated  for  a terminated 
lying  on  the  ocean  floor. 


See  Figure  6.2  for  geometry  employed 


4.  The  Two-Half-Space  Model  and  Lateral-Wave  Transmission 
4.1  Introduction 

In  their  modal  analysis  of  electromagnetic  wave  propagation  in  a litho- 
sphere composed  of  granitic  rock  (j=  10  ° Si/m)  and  bounded  above  by  sea 
water  (o  ■ 4 Si/m)  and  below  by  the  mantle  (o  = « ) Frieman  and  Kroll  [1973] 
concluded  that  with  horizontal-electric-dipole  excitation  and  reception  at 
points  on  the  ocean  floor,  the  received  field  was  very  little  affected  by 
the  thickness  of  the  lithosphere  (d„  "v  10  km)  or  by  the  depth  of  the  water 

(d  ^ 5 km)  in  the  frequency  range  from  0.1  to  10  Hz.  This  means  that  the 

w 

electric  field  at  the  point  of  reception  is  governed  by  a single  dominant 
mode  which  corresponds  to  the  solution  for  transmission  along  the  boundary 
between  two  half-spaces,  the  ocean  and  the  lithosphere.  This  conclusion 
is  confirmed  by  the  careful  interpretation  by  Bubenik  and  Fraser-Smith  [1978] 
of  their  analysis  of  the  field  of  a vertical  magnetic  dipole  located  in 
sea  water  midway  between  the  air  above  and  the  sea  floor  (lithosphere)  below. 

Bubenik  and  Fraser-Smith  compared  the  contributions  to  the  total  field  by 
direct  transmission  through  the  sea  as  if  this  were  infinite  in  all  directions 
with  those  by  air-surface  and  sea-floor  lateral-wave  modes.  They  concluded 
that  for  a dipole  near  the  air  surface,  the  air-surface  mode  dominates;  for 
a dipole  near  the  sea  floor,  the  sea-floor  mode  dominates.  Since  a vertical 
magnetic  dipole  is  equivalent  to  four  horizontal  electric  dipoles  arranged 
in  a square,  its  field  is  a superposition  of  the  fields  of  the  four  dipoles;  , 
each  like  the  field  of  a single  horizontal  electric  dipole.  The  superimposed  fields 
of  the  very  closely  spaced  dipoles  lead  to  a large  reduction  in  the  resultant  field 
due  to  cancellation  but  this  in  no  way  alters  the  nature  of  the  field  with 
respect  to  its  modal  structure.  This  means  that  if  the  sea-floor  lateral-wave 
mode  dominates  for  a vertical  magnetic  dipole  in  sea  water  near  the  sea  floor, 
the  same  is  true  for  a horizontal  electric  dipole  at  the  same  location. 
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Extensive  numerical  calculations  by  King  and  Sandler  [1977]  based 

on  the  exact  general  theory  of  Banos  (1966]  have  shown  that  transmission 

between  horizontal  electric  dipoles  in  a dissipative  half-space  near  its 

boundary  with  air  is  predominantly  by  the  so-called  lateral  waves.  This  is 

true  for  a dissipative  half-space  that  has  the  properties  of  sea  water, 

9 

lake  water,  or  dry  earth  over  a frequency  range  from  1 to  10  Hz.  For 
radial  distances  from  o = 1 m to  100  km  and  more,  the  dominant  component 
of  the  electric  field  at  a depth  z in  the  medium  due  to  a transmitting  dipole 

at  the  depth  d is  the  radial  one,  i.e.,  E . While  there  are  short  ranges 

p 

where  the  component  E,  is  somewhat  greater  than  E , this  difference  is 

9 P 

never  large.  On  the  other  hand,  there  are  large  ranges  including  especially 

the  greater  distances  where  E greatly  exceeds  E . It  follows  that  for 

P $ 

practical  transmission  over  large  distances  E is  the  only  satisfactory 

P 

component.  It  has  also  been  shown  by  direct  comparison  of  the  numerical 
calculations  from  the  exact  integrals  with  a set  of  approximate  formulas 
due  to  Banos  [1966]  that  these  latter  are  roughly  correct  at  least  in  order 
of  magnitude  far  beyond  the  ranges  of  validity  specified  by  Banos.  Specifically, 
the  graphs  of  the  electric  field  as  a function  of  distance  at  a given  frequency 
as  calculated  from  3anos'  near-,  intermediate-,  and  far-field  formulas  can  be 
extended  to  points  of  intersection.  The  sequential  combination  of  the  three 

graphs  provides  a continuous  curve  of  I E 1 against  distance  from  the  source  3 

‘ P ' 

(but  with  discontinuities  in  slope  at  the  intersections  of  the  near-  and  inter- 
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mediate-field  graphs  and  of  the  intermediate-  and  far-field  graphs)  that 
follows  fairly  closely  the  smooth  graph  computed  from  the  exact  integrals. 

This  general  agreement  is  very  significant  because  it  justifies  the  use  of 
Banos'  simple  formulas  to  gain  an  understanding  of  the  physical  nature  of 
the  waves  propagating  from  the  source  to  the  receiver  and  of  the  significance 
of  the  numerous  parameters  and  variables  that  combine  to  determine  the  field 
at  the  point  of  observation. 

4.2  Banos'  Approximate  Formulas  and  Their  Interpretation 

The  coordinates  and  parameters  used  in  the  representation  of  the  field 

in  a region  consisting  of  a half-space  of  sea  water  (Region  1,  z > 0)  and  a 

half-space  of  the  lithosphere  (Region  2,  z < 0)  are  shown  in  Fig.  4.1  . The 

wave  number  in  sea  water  with  the  time  dependence  e lwt  is* 

2 1/2  1/2 

kl  " S1  + ictl  = (“  u0£l  + ^Vl5  4 (4.1) 

The  wave  number  for  the  lithosphere  is 


k-  = 3„  + ia  „ = (a  + iaiunJ0)' 


(4.2) 


As  indicated,  the  sea  water  behaves  like  a good  conductor  over  all  useful 
frequencies.  The  same  is  true  for  the  lithosphere  at  low  frequencies.  At 
high  frequencies  the  lithosphere  may  behave  like  a dielectric  so  that  the 
more  general  formula  for  the  wave  number  is  retained. 

Bafios'  approximate  formulas  for  the  electric  field  of  a horizontal  electric 
dipole  at  the  location  shown  in  Fig.  4.1a  are  given  by  (7.10),  (7.32)  and  (7.45) 
in  his  book  [Banos,  1966].  With  small  changes  in  the  notation  and  the  substi- 


tution of  k“/iuiUQ  for  a ^ , the  formulas  for  the  radial  component  of  the  electric 


field  are: 


Near  Field: 


iu)Un  „ j,  i [k.p  + k (d+z)] 
U cos  9 2 1 

2 3 e 

2rk^  o 


(4.3) 


Intermediate  Field: 


*The  form  on  the  right  in  (4.1)  assumes 


1 for  sea  water. 


Fig.  4.1a  Two-half-space  Model  for  Lateral-Wave  Transmission 
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io)Ur.k^  cos  5 k.  . i [k?p  + k.(d+z)] 

E.  A 7-=- [1  + t — (iik  p/2) 1/ Z ] e " 1 (4.4) 

?n\c~  0 k. 


Far  Field: 


- • MU0  cos  > Hk2p  + kx(d+z)] 

Llp  ’ 2wk.  2 e 

2 0 


(4.5) 


Tae  formulas  for  the  other  components  are  similar,  but  since  is  the  only 
generally  useful  one,  it  is  sufficient  to  consider  it.  The  ranges  in  which 
each  of  the  three  formulas  is  a good  approximation  have  been  specified  by 
Banos.  For  sea  water  they  are  extensive,  but  even  for  this  they  do  not  actuallv 
join.  Numerical  calculations  from  the  exact  integrals  have  shown  the  three 


formulas  may  be  used  successively  as  the  radial  distance  is  increased  with 
the  transition  from  the  one  to  the  other  occurring  where  the  magnitudes  of  the 
field  are  equal.  At  and  near  these  transitions  from  the  one  formula  to  the 
other  the  approximations  are  quantitatively  inadequate  but  qualitatively 


correct . 


Tlie  lateral-wave  nature  of  the  field  is  evident  in  all  three  formulas. 

Note  that  they  all  have  the  same  exponential  term  which  is  readily  interpretted . 
It  indicates  a traveling  wave  originating  at  the  source  dipole,  (p=0,  p=0,z=d), 
traveling  downward  in  the  sea  water  to  the  boundary  surface  , then  traveling 
radially  outward  in  the  lithosphere  a distance  p , and  finally  traveling  upward 
from  the  surface  to  the  receiver  at  the  point  of  observation  at  P(p,5,z).  This 
involves  travel  in  the  sea  water  (Region  1)  over  a total  distance  z + d and 
travel  in  the  lithosphere  (Region  2)  over  the  distance  o rFig.  4,1b).  In  the  sea 
the  attenuation  and  phase  constants  are,  respectively,  a1  and  3^;  in  the  litho- 
sphere they  are  and  If  the  real  and  imaginary  parts  of  the  exponents 

are  separated,  the  common  exponential  term  has  the  expanded  form 


-[c^P  + a^(d+z)]  i[30p  + 3^(d+z) 

e e 


(4.6) 
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Clearly  there  is  attenuation  in  both  regions  but  since  the  conductivity  of  the 
lithosphere  (Region  2)  is  much  smaller  than  that  of  the  salt  water  (Region  1) 
it  follows  that  <<  a^.  For  a given  range  p an  important  quantity  in  deter- 
mining the  amplitude  is  the  total  vertical  distance  z + d.  The  closer  the  trans- 
mitting and  receiving  antennas  are  to  the  boundary,  the  larger  will  be  the  ampli- 
tude of  the  field.  However,  if  a^(d+z)  < 1,  g^^W+z)  ^ ^ even  for  large 

values  of  d + z.  This  can  be  accomplished  by  choosing  a low  enough  frequency  so 
that  is  small  even  though  is  as  large  as  4 S/m.  Insofar  as  the  exponential 

term  is  concerned,  small  distances  z and  d are  desirable  as  well  as  a frequency 
that  is  not  too  high.  For  a specified  range  the  contribution  to  the  attenuation 

“Ct  0 

by  the  factor  e 2 can  be  reduced  only  by  making  ct0  smaller.  With  a given  conduc- 

tivity  of  the  lithesphere  a„  can  be  decreased  only  by  lowering  the  frequency. 

If  the  exponential  term  does  not  dominate, a significant  field  can  be 

maintained  at  great  distances.  In  order  to  understand  its  dependence  on  the 

distance  and  frequency  it  is  necessary  to  examine  the  other  factors  that 

determine  the  amplitude.  The  radial  distance  appears  in  all  three  formulas 

(4.3)  - (4.5)  but  in  a different  manner  in  each.  In  the  near-field  formula 

3 

the  amplitude  decreases  as  l/o  , a prohibitively  rapid  rate  of  decrease. 

In  the  far-fiela  formula  the  rate  of  decrease  is  as  l/o~  , a great  improvement 

over  the  near  field.  Finally,  in  the  intermediate  range  the  rate  of  decrease 

with  radial  distance  varies  between  l/o  and  l/»c  - which  is  by  far  the  most 

favorable.  The  comparable  dimensionless  quantities  for  the  three  ranges  are: 

-3  2 3 2 

(k^o)  , k“/k^c  , and  l/k^k7o  in  the  order  near  field,  intermediate  field,  and 

far  field.  With  a proper  choice  of  frequency  for  the  parameters  of  the  media 
involved,  the  largest  magnitude  of  the  radial  electric  field  should  be  obtained 
'-•hen  o is  in  the  intermediate  zone  near  its  boundary  with  the  far  zone.  This 
clearly  excludes  very  low  frequencies  for  which  the  near-field  formula  is  appli- 
cable for  the  relevant  radial  range. 

A better  understanding  of  the  complete  picture  including  the  exponential 
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and  inverse  distance  dependences  can  be  obtained  from  Fig.  4.2  which  displays 

contours  of  constant  IE  (in  d3  referred  to  1 volt  per  meter)  as  a function  of 

P 

the  radial  distance  and  the  frequency  for  a horizontal  electric  dipole  at  a 
depth  a = 15  cm  in  water.  The  point  of  observation  is  at  the  depth  z = 30  cm  in 
the  water.  The  half-space  bounding  the  water  is  air  filled.  In  the  previous 
notation,  the  water  is  Region  1,  the  air  is  Region  2.  In  this  case  012  = 0 and 
$2  = 2tt / \ 2 where  X2  is  the  wavelength  in  air.  The  real  effective  permittivity  of 
the  water  is  e , = 30;  the  real  effective  conductivity-  in  Siemens  per  meter  is 
a ^ = 0.004,  0.04,  0.4  and  4 in  the  diagrams  from  left  to  right.  The  diagram 
on  the  left  applies  when  the  antennas  are  in  lake  water  of  high  purity,  the  dia- 
gram second  from  the  left  is  for  polluted  lake  water,  the  diagram  on  the  extreme 
right  for  sea  water,  The  contours  were  obtained  from  numerical  computations 
of  the  electric  field  using  the  exact  integrals. 

In  the  diagram  on  the  left  the  conductivity  is  sufficiently  low  so  that 

the  exponential  attenuation,  e is  ll0C  dominant  at  any  frequency  with 

z + d = 45  cm.  In  terms  of  the  approximate  formulas  of  Banos,  the  right-hand 

third  of  the  diagram  is  for  low  frequencies  and  the  near-field  formula  with 

,3 

its  l/o  amplitude  factor,  the  left-hand  third  is  for  high  frequencies  and  the 
far-field  formula  with  its  1/p ^ amplitude  factor.  The  transition  between  them 
at  intermediate  frequencies  is  for  the  intermediate-field  formula.  Note  that  at 
any  fixed  distance,  for  example  c = 8 km,  is  very  much  greater  at  the 

high  frequencies  than  at  the  low  ones.  Specif ically , at  the  high  frequencies 
E,  = -120  dB,  in  the  low-frequency  range  E.  = -200  dB.  In  the  second 

-*-C  lo 

diagram  from  the  left  the  conductivity  of  the  water  is  ten  times  greater  but 
still  small  enough  to  make  the  exponential  attenuation  unimportant.  The  contours 
of  constant  electric  field  are  quite  similar  to  those  on  the  left  with  only  a 
small  decrease  in  amplitude  in  the  high-frequency  range,  a somewhat  larger  one 
in  the  low-frequency  range.  In  the  third  diagram  from  the  left  the  conductivity 
is  multiplied  by  another  factor  ten  and  this  makes  it  large  enough  so  that  the 
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exponential  attenuation  begins  to  play  a significant  role,  but  only  at  the 
nigner  frequencies.  The  rapid  decrease  in  amplitude  at  frequencies  above  10^  Hz 
is  due  to  the  exponential  factor  e 3-^2+c^^  with  z+d  = 45  cm.  Note,  however, 
that  even  with  the  exponential  attenuation  significant,  the  field  in  the  high 
frequency  range  is  still  much  greater  at  any  fixed  distance  than  in  the  low 
frequency  range.  The  optimum  frequency  is  clearly  in  the  vicinity  of  10/  Hz 
where  the  contours  show  a maximum  distance  for  any  specified  value  of  E j. 

In  Che  graph  on  the  right  for  sea  water  the  exponential  attenuation  has 
become  the  major  factor  at  the  higher  frequencies.  The  radial  electric 

field  decreases  to  extremely  small  values  at  all  distances  when  the  frequency 
g 

is  increased  to  10  Hz  or  higher.  However,  at  any  fixed  distance  the  field 
has  a maximum  near  f = 10 / Hz  and  this  maximum  is  significantly  greater  than 
the  field  at  the  same  distance  at  any  of  the  lower  frequencies.  Thus,  when 
z + d = 45  cm  the  field  at  any  radial  distance  p is  greater  at  a frequency 
between  5 x 10^  and  5 x 10 7 than  at  any  frequency  below  f = 10^  Hz. 

If  the  total  distance  of  travel  in  the  sea  water  is  increased  from  z+d 
= 45cm  the  maximum  field  at  any  fixed  radial  distance  moves  toward  lower 
frequencies.  When  this  reaches  about  103  Hz  the  exponential  attenuation  will 
effectively  eliminate  the  maxima  in  the  high  frequency  range  and  the  field 
in  the  low-frequency  range  will  be  larger  than  that  anywhere  in  the  high-fre- 
quency range.  The  choice  of  frequency  for  lateral  wave  transmission  is  seen 

mmrX  ( z+d ) 

to  depend  critically  on  the  quantity  i^(z  + d)  and  the  exponential  e 1 
linen  the  region  of  low  conductivity  is  not  air  with  a,  = 0 but  the  lithosphere 

-8  -4 

with  o0  = 10  to  10  .the  additional  exponential  attenuation  by  the  factor 
e must  also  be  considered.  This  requires  a complete  recalculation. 

However,  the  basic  properties  of  wave  propagation  with  horizontal-dipole 
excitation  near  an  interface  are  best  understood  from  the  simpler  case  with 
a perfect  dielectric  like  air  as  the  half-space  with  small  conductivity. 

3efore  proceeding  to  a study  of  lateral  ’.rave  propagation  along  the 
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surface  of  two  half-spaces  with  the  properties  of  sea  water  and  the  lithosphere, 
it  is  well  to  point  out  that  the  simple  picture  illustrated  in  Fig.  4.1b  of 
lateral  waves  traveling  vertically  downward  from  the  source  dipole  a distance 
d to  the  boundary  surface,  along  the  surface  a radial  distance  p to  the  point 
below  the  receiver,  and  then  vertically  upward  a distance  z to  the  receiving 
dipole,  is  an  idealization  that  is  an  excellent  approximation  with  a highly 
conducting  medium  like  sea  water.  With  poorer  conductors  the  optimum  angle 
downward  from  the  source  to  the  surface  is  the  critical  angle  at  which  the 
incident  wave  is  neither  transmitted  nor  totally  reflected  but  travels  along 
the  surface.  With  sea  water  this  angle  is  very  close  to  90°;  for  media  with 
lower  conductivities  and  permittivities  the  angle  can  be  considerably  smaller 
than  90°.  A discussion  and  determination  of  this  angle  is  given  by  Staiman 
and  Tamir  [1966].  For  a maximum  transmitted  signal  the  dipole  should  not 
be  parallel  to  the  interface  but  perpendicular  to  the  direction  of  the  critical 
angle.  With  sea  water  the  difference  between  the  two  orientations  is  negligible. 
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4.3  Lateral -Wave  Transmission  Along  the  Lithosphere-Ocean  Boundary 

Lateral-wave  transmission  along  the  boundary  between  air  and  water  involves 
exponential  attenuation  only  along  the  two  vertical  parts  of  the  path,  both  in 
the  salt  water.  This  has  the  form  e-0*1 (z+^)  where  j the  depth  of  the  horizon- 
tal dipole  in  the  sea  water  and  z is  the  depth  of  the  point  of  observation.  When 
the  transmission  is  along  the  boundary  in  air  it  is  loss  free  since  e^20  = e*""-" 
with  k,,  = + ia0  and  = 0.  When  Che  half-space  of  air  is  replaced  by  a naif- 

space  of  rock  with  the  real  effective  permittivity  eg?  and  the  real  effective 

conductivitv  a the  wave  number  k„  = 8„  + ict„  is  complex  so  that  e^20  = 
e2  2 2 2 

ei32P  • Thus  there  is  exponential  attenuation  along  the  main  path  of 
propagation  in  the  lithosphere  in  addition  to  the  decrease  in  amplitude  governed 
by  the  complicated  functional  dependence  on  the  inverse  powers  of  the  radial 
distance  p. 

Computations  have  been  made  from  the  exact  integrals  for  lateral-wave 
transmission  from  a horizontal  electric  dipole  with  unit  moment,  Id£  = 1. 

The  dipole  is  at  the  depth  d in  the  sea  water  and  the  field  is  calculated  at 
the  radial  distance  p at  the  depth  z = d.  (Note  that  in  the  computations  for 
air  and  sea  water,  z = 2d  = 0.3m)  . The  material  parameters  used  in  the  cal- 
culations of  the  electric  field  are  e , = 80e„  and  o , = 4 Si/m  for  Region  1, 

el  0 el 
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the  sea.  For  the  lithosphere,  Region  2,  t 7 = 16eq  and  ranges  from  4 x 10 

Si/ra  to  4 x 10  Si/m.  Initially  z = d = 0.15m  . Computations  of  !E  ' (which 

1 0 

is  the  dominant  component  in  lateral-wave  transmission)  have  been  made  over 

9 

radial  distances  from  p=  0.1  to  p = 100km  at  frequencies  from  10  to  10  Hz. 

Graphs  showing  1 E ! as  a function  of  the  radial  distance  with  the  frequency  as 

P 
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the  parameter  are  shown  in  Fig.  4.3a  for  o = 4 x 10  and  4 x 10  Si/m  and 
in  Fig.  4.3b  for  a , * 4 x 10  ^ 4 x 10  \ 4 x 10  \ and  4 x 10  ^ Si/m.  From 

these  graphs  the  contours  of  constant  |E  ! shown  in  Figs.  4.4a  - 4.4c  have  been 

P 

constructed.  These  show  E ! as  a function  of  both  the  radial  distance  and 

P 

the  frequency  in  the  manner  carried  out  in  the  comparable  Fig.  4.2  for  air  and 


water . 


-4.17- 


The  set  of  contours  on  the  left  in  Fig.  4.4a  is  for  sea  water  and  a litho- 

—8 

sphere  with  the  low  conductivity  a = 4 x 10  Si/m.  It  resembles  the  correspon- 
ding set  of  contours  on  the  extreme  right  in  Fig.  4.2  for  sea  water  and  air  with 
the  conductivity  c?e?  = 0.  In  both  sets  of  graphs  the  exponential  attenuation 
e -M  dominates  the  high-frequency  quarter  on  the  left.  Here  the  graphs  show 

a rapidly  decreasing  radial  range  for  any  given  value  of  j E | as  the  frequency 

P 

is  increased.  The  rate  of  decrease  for  sea  water  and  air  is  even  more  rapid  than 

that  for  sea  water  and  the  lithosphere  because  (z+d)  = 0.45m  with  the  former,  0.3m 

with  the  latter.  On  the  other  hand,  attenuation  by  e 1 is  insignificant  in 

the  low-frequency  third  of  the  contours  on  the  right  in  each  diagram.  This  is 

because  at  low  frequencies  <*2  is  very  small.  For  air,  a2  = 0,  e-020  = 1,  and 

the  contours  for  fixed  values  of  are  virtually  horizontal  in  the  low-frequency 

third  of  the  right-hand  set  of  contours  in  Fig.  4.2.  This  indicates  that  the 
radial  range  for  a specified  field  is  almost  independent  of  the  frequency  in  the 
low-frequency  range.  For  the  lithosphere,  U2  ^ 0 and  e a“°  varies  significantly 

with  both  the  frequency  and  the  radial  distance.  As  a consequence,  the  contours 
for  given  values  of  E at  low  frequencies  range  from  horizontal  lines  when  c ^ 
is  small  to  steeply  increasing  graphs  when  a ~ is  large.  The  increase  in  radial 
range  as  the  frequency  is  reduced  and  the  conductivity  increased  is  clearly  shown 
in  the  sequence  of  graphs  in  Figs.  4.4b  and  4.4c. 

Perhaps  the  most  significant  characteristic  of  the  contours  of  constant 

, E ; for  sea  water  and  air  (on  the  extreme  right  in  Fig.  4.2)  and  for  sea  water 

P 

and  the  lithosphere  (Figs.  4.4a  - 4.4c)  is  the  relative  maximum  in  the  radial 

range  which  occurs  in  all  diagrams  at  frequencies  between  f = 1 and  10  MHz  so 

_3 

long  as  a ,,<  4 x 10  Si/m.  Particularly  notworthy  is  the  fact  that  the  rela- 
tive maximum  value  of  p for  a given  E_  is  much  greater  than  the  corresponding 

distance  for  the  same  value  of  E : at  very  low  frequencies  when  a ^ * 4 x 10 

“8 

Si/m.  This  can  be  seen  from  Figs.  4.4a  and  4.^b  for  both  a ^ = 4 x 10 

t , « 4 x 10  Si/m.  For  example,  when  p , 5 4 x 10  Si/m,  the  value 


and 
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E ; * -175  dB  occurs  at  c 3 40  km  when  the  frequency  is  between  10^  and  10^  Hz 

but  o < 0.3  km  for  all  frequencies  from  10  to  10"*  Hz.  Similarly  [E  [ = -225  dB 

occurs  at  o near  500  km  when  the  frequency  is  between  10^  and  10^  Hz  but  at 

4 -7 

2 km  for  all  frequencies  between  10  and  10  Hz.  When  a „ = 4 x 10  Si/m  the 

e2 
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graphs  are  very  similar  to  those  with  a 0 = 4 x 10  Si/m  with  a reduction  to 
p ■ 100  km  in  the  radial  range  at  the  maximum  between  f = lO*3  and  10 ' Hz  and 
virtually  no  change  in  the  short  range  at  low  frequencies.  However,  when  og0 
is  as  large  as  4 x 10  b Si/m  the  maxima  in  the  radial  range  near  f = 10^  Hz  are 

substantially  depressed  but  the  radial  distance  for  the  same  values  of  Je  ' at 

P 

very  low  frequencies  remains  smaller  so  long  as  p < 75  km.  Thus,  | E I = -225 

o 

dB  occurs  at  p * 33  km  at  the  maximum  with  the  frequency  between  10°  and  10 ‘ Hz 

3 

and  drops  to  p = 2 km  at  f = 10  to  10  Hz.  For  radial  distances  greater  than 

about  100  km,  a given  value  of  |E  | is  found  at  substantially  greater  distances 

P 

at  the  very  low  frequencies  than  at  the  relative  maximum  between  f = 10°  and  lo" 

Hz.  This  is  a consequence  of  the  radial  exponential  attenuation  factor  e 

P 7 

which  is  small  enough  at  frequencies  as  high  as  10  and  10  Hz  to  depress  the 
relative  maximum  more  and  more  as  p is  increased.  This  effect  is  dramatically 

shown  in  the  four  sets  of  graphs  in  Figs.  4.4b  and  4.4c  in  which  the  conductivity 

”6  * 3 

of  the  lithosphere  increases  from  4 x 10  on  the  left  in  Fig.  4.4b  to  4 x 10 

n 

Si/m  on  the  right  in  Fig.  4.4c.  It  is  seen  that  E \ at  the  maximum  near  f = 10 

p 

Hz  is  reduced  more  and  more  while  |E_J  at  very  low  frequencies  and  sufficiently 
great  radial  distances  is  increased. 

The  contour  diagrams  in  Figs.  4.4a,b  lead  to  the  following  general  conclu- 
sions regarding  lateral-wave  transmission  along  a plane  boundary  between  sea  water 
and  the  lithosphere  when  the  source  is  a horizontal  electric  dipole  in  the  sea 
at  the  small  distance  d = 15  cm  from  the  bounding  surface  and  the  point  of  obser- 
vation is  at  the  same  distance  z = d from  the  surface  at  a radial  distance  o : 
a)  If  the  conductivity  of  the  lithosphere  is  greater  than  * 4 x 10  ° Si/m, 


transmission  over  distances  greater  than  50  km  experiences  the  lowest  possible 
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attenuacion  at  frequencies  below  1 kHz.  This  is  true  for  any  values  of  z = d 
>15  cm. 

b)  If  the  conductivity  of  the  lithosphere  is  significantly  smaller  than  4 x 10  ^ 
Si/m,  there  is  an  optimum  frequency  between  1 and  10  MHz  at  which  the  radial 
range  has  a relative  maximum  that  may  be  substantially  greater  than  the  radial 
distance  for  any  value  of  |E  | at  other  frequencies  including  the  very  low. 
Furthermore,  an  effective  directional  array  with  significant  gain  can  be  con- 
structed for  both  the  transmitter  and  the  receiver  for  use  at  frequencies  in 
the  1 to  10  MHz  range. 

c)  The  conclusion  under  b)  assumes  the  transmitting  and  receiving  antennas  to 

be  in  the  sea  water  but  very  near  the  boundary  surface  with  the  lithosphere 

(d  = z = 0.15m).  It  is  to  be  expected  that  as  the  distances  d and  z from  the 

antennas  to  the  surface  is  increased,  the  exponential  attenuation  e-'M(z+t^ 

will  become  dominant  at  lower  and  lower  frequencies.  This  must  simultaneusly 

depress  the  relative  maximum  which  occurs  between  f = 1 and  10  MHz  when  (z+d)  = 

30  cm  and  move  it  to  lower  frequencies.  In  order  to  investigate  this  effect 

quantitatively .computations  were  made  of  ' with  z = d = 1.5  m over  the  same 

range  of  conductivities  for  the  lithosphere  as  used  in  Figs.  4.3a  and  4.3b. 

These  results  are  shown  in  Figs.  4.5a  and  4.5b.  From  them  the  contour  diagrams 

in  Figs.  4.t>a  - 4.6c  were  constructed.  They  correspond  to  the  diagrams  in 

Figs.  4.4a  - 4.4c.  It  is  seen  that  when  a , < 4 x 10  3 Si/m  the  relative 

maxima  in  the  radial  range  have  moved  from  near  f = 10^  Hz  wit1'  z * d = 0.15m 

to  near  f = 103  Hz  with  z = d = 1.5m  and  that  they  are  substantially  smaller. 

—8 

For  example,  with  a = 4 x 10  Si/m  the  maximum  range  for  |E  » -225  dB  is 
p = 10  km  with  z = d = 1.5m  instead  of  p v 500  km  with  z = d = 0.15m.  There  is 
no  corresponding  decrease  in  the  range  at  low  frequencies  since  there  the  expo- 
nential attenuation  is  small.  It  follows  that  unless  both  the  transmitting  and 
receiving  antennas  are  very  close  to  the  ocean-lithosphere  boundary,  the  greatest 
radial  range  is  achieved  at  low  frequencies. 
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Magnitude  of  E p in  Sea  Water  Bounded  by  the  Lithosphere 
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d)  There  is  a potentially  useful  low-attenuation  window  in  the  frequency  range 
between  0.1  and  10  MHz.  It  exists  only  when  the  conductivity  of  the  lithosphere 
is  sufficiently  low,  < 4 x 10  ° Si/m.  Furthermore,  it  is  available  only 

when  the  transmitting  and  receiving  antennas  are  very  close  to  the  interface 
between  the  ocean  and  the  lithosphere.  When  (z-M)  = 0.3m  the  maximum  radial 

range  occurs  when  the  frequency  is  near  10^  Hz.  As  the  distance  from  the 
boundary  is  increased  the  frequency  for  the  maximum  range  decreases.  It  is  near 
f = 10^  Hz  when  (z+d)  = 3 m.  The  magnitude  of  the  maximum  range  also  decreases 
i as  (z+d)  is  increased. 
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5.  LABORATORY  MODELLING  AND  ELECTRODYNAMIC  SIMILITUDE 

5.1  Introduction 

Laboratory  modelling  has  been  used  to  simulate  or  demonstrate  the 
behavior  of  large  scale  electromagnetic  systems  by  scaling  dimensions  and 
frequencies  while  keeping  the  same  shape  of  the  model  and  the  system  it  is 
to  simulate.  For  example,  a full  scale  linear  antenna  of  length  Lg  operating 
in  air  at  a frequency  f can  be  simulated  by  measurements  on  a similar  shaped 
antenna  of  length  = Le/n  provided  the  model  is  operated  at  a frequency 
fm  = n fe  in  air.  For  more  general  lossy  media,  a knowledge  of  the  prin- 
ciples of  electrodynamic  similitude  and  limitations  is  needed  for  the  use 
of  models  and  interpretation  of  data.  We  examine  relations  for  application 
to  some  propagation  problems. 

Millington  (1948,  1949a, b)  discovered  the  "recovery  effect"  on  the 
amplitude  of  LF  and  VLF  ground  waves  as  they  traversed  a land-sea  boundary. 

A simple  model  was  used  to  demonstrate  the  effect,  wood  to  represent  the 
land  and  copper  sheeting  to  represent  the  sea.  A small  microwave  source  was 
moved  on  a small  cart  away  from  the  source  on  wood  onto  the  metal  whence  the 
signal  amplitude  increased  markedly  at  the  "coast-line,"  (recovery)  and  then 
attenuated  slowly  along  the  copper  "sea".  Quasi-empirical  relations  were 
developed  in  his  publications.*  Pressey  (1956)  developed  quasi-empirical 
relations  for  the  phase  shift.  The  matter  was  studied  theoretically  and 
experimentally  by  King  and  Maley  (1966)  and  King,  Maley,  and  Wait  (1967) . 

*A  group  of  us  in  February,  1950,  witnessed  Millington's  demonstration 


at  Marconi's  arranged  through  the  courtesy  of  the  late  Sir  Robert  Watson  Watt. 
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The  model  was  a tank  filled  with  water  to  represent  land  and  aluminum  sheets 
to  represent  the  ocean;  data  were  obtained  at  4 GHz.  While  scaling  was  not 
exact,  the  measured  electric  field  strength  was  reasonably  well  in  accord 
with  prediction. 

Model  studies  of  propagation  in  an  idealized  crustal  waveguide  were 
made  experimentally  and  theoretically  by  Brown  and  Gangi  (1962,  1963)  and 
Gangi  (1966) . We  shall  discuss  their  results,  arranged  to  demonstrate 
similitude,  a bit  later,  for  a uniform  slightly  lossy  guide  represented  by 
saline  water. 

Iizuka  (1967,  1968,  1969)  studied  the  impedance  of  an  antenna  and  the 
field  variations  of  waves  propagated  in  a lithospheric  waveguide  for  a guide 
having  variable  conductivity  with  depth.  A gel  of  agar-agar  with  diffusion 
of  sodium  chloride  was  arranged  to  provide  the  desired  profile.  Similitude 
in  size,  frequency,  etc.,  was  accomplished  and  results  will  be  discussed  a 
bit  later,  also.  Others  have  developed  the  similitude  relations. 

5.2  Principles  of  electrodyr.amic  similitude 

The  principles  of  electrodynamic  similitude  were  developed  by  Stratton 
(1941)  and  more  completely  delineated  by  Sinclair  (1948) . Brown  and  Gangi 
(1962,  1963)  and  Gangi  (1966)  gave  the  resulting  equations,  based  on 
Stratton  (1941)  relating  the  scaling  of  frequency,  sizes,  and  electrical 
parameters  of  the  model  system  in  terms  of  those  quantities  of  the  original 
system.  Keller  and  Fritschkneckt  (1966)  gave  a brief  discussion  for 
scaling  media  of  large  loss  tangent.  Iizuka  (1968)  developed  the  relations 
while  employing  generalized  orthogonal  coordinates  for  the  identities  of 


Maxwell's  equations  in  the  original  and  model  systems.  Fritschkneckt  (1971), 
basing  his  work  in  part  on  that  of  Sinclair  (1948) , used  rectangular  coordin- 
ates as  did  Sinclair  (1948). 


Since  we  are  dealing  with  systems  immersed  in  the  earth,  such  as  the 
lithosphere,  we  use  the  subscript  "e"  on  related  quantities;  the  subscript 
"m"  will  denote  those  quantities  for  the  model. 


The  original  and  model  media  are  characterized  as  isotropic,  homogeneous, 

or  inhomogeneous,  with  electrical  constants  of  permittivity  c = e0  er  (where 
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eQ  is  approximately  10  /36n  F/m  and  er  is  the  relative  dielectric  constant, 

a permability  y (which  we  will  ultimately  assume  to  be  that  of  free  space 
uD  = 4it  x 10-7  H/m)  and  a conductivity  a in  Si/m.  The  unit  of  length  is  r. 

In  order  that  Maxwell's  equations  remain  invariant  from  the  original  to  model 
system,  we  shall  scale  the  electric  field  strength  E,  the  magnetic  field 
strength  H,  and  time.  Thus,  for  the  fields: 


% = Ee  (5.1) 

= qH  He  <5-2> 

By  showing  the  invariance  of  Maxwell's  equations,  we  arrive  at  the  rela- 
tions for  Um,  em  and  dm. 


Um  q. 


re 


e = 
m 


fe 

^E  rm 

^E  rm 


where  T is  the  unit  of  time  such  that  Tn\/Te 


(5.3) 

(5.4) 

(5.5) 

fe/fm  = “e^m  where  f is  the 
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wave  frequency  and  oo  the  angular  frequency.  The  equations  (5.1  through  5.3) 
are  cast  in  the  form  of  Iizuka  (1968)*;  as  mentioned  by  Fritschkneckt  (1971), 
they  are  termed  by  Sinclair  (1948)  as  "absolute"  models,  for  which  qE  and  qH 
are  known  separately  and  not  merely  their  ratio.  If  a specific  value  is  not 
assigned  to  qE  or  q^,  the  model  is  termed  "geometric"  by  Sinclair  (1948) . 


Multiplying  equation  (5.3)  by  (5.4),  one  obtains: 


r_  w_  /ume_  = r oj  /u  £ 
m m mm  ee  e e 


(5. 6. a) 


k_'r  = k ' r 
mm  e e 


(5.6.b) 


r /\  ' = r /X  ' 
nr  m e e 


(5.6.c) 


where  k ' = wj/ue  = 2ir/X  ' and  k ' = w </iT~e7  ~ 2ir/X^ ' ; i.e.  , the  k ' 
m nr  m m m e eee  e m 


and  k ' are  the  wave  numbers  in  loss-less  dielectrics.** 
e 


From  (5.4)  and  (5.5): 


(5. 7. a) 


U)  € OJ  € 

mm  e e 


60 o X 60*7  X 

o = o _ — _ e o 

m e e £ 

rm  re 


(5.7.b) 


where  d is  also  known  as  the  loss-tangent  tan  5 and  X„  is  the  free 
m,e  m,e° 


space  wavelength  (meters). 


'i.e. , equations  (IV-13)  through  (TV-15)  of  Iizuka  (1968) 


♦♦Equations  (5.6)  through  (5.9)  correspond  to  Iizuka 's  (1968)  equations 
(TV-16)  through  (IV-19) . 
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Again,  from  equation  (5.3) 


— r a u = r w u 
in  m m e e e 


(5.8) 


^E 


which  can  be  changed,*  by  using  (5.6) , to  show: 


?E  fee  l^m 


q„  e ' U 1 £ > U 

me  rm  e 


(5.9) 


This  relation  is  used  when  showing  the  scaled  impedance  relation  of  original 
and  modelled  antennas. 

lizuka  (1968)  preferred  as  convenient  to  express  the  similitude  condi- 
tions as  follows:  equation  (5.6)  states  the  identity  of  electrical  sizes  in 
terms  of  the  loss-less  wavelength  between  original  ("e")  and  model  ("m") 
systems,  and  equation  (5.7)  states  the  identity  of  loss-tangents  cf  the  two 
systems. 

Brown  and  Gangi  (1962,  1963)  and  Gangi  (1966)  gave  relations  similar  to 
those  above,  viz.: 


f u e r -f  u e r 
m m m m e e e e 


fm  um  °m  V = fe  ue  ee  V 


(5.10) 


(5.11) 


If  one  takes  the  ratio  of  (5.11)  to  (5.10),  the  identity  of  loss-tangents 
in  equations  (5.7)  result. 


•Equations  (5.6)  through  (5.9)  correspond  to  lizuka' s (1968)  equations 
(IV-16)  through  (IV-19) . 


5.6 


i 
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With  the  assumption  that  U_  = U = 


Vo,  they  go  to  show: 


a f 
m _ rm 


,/Srm  _ De  /fra i 


(5.12) 


£m 

fd 


(5.13) 


where  De  is  the  depth  of  the  original  lithospheric  waveguide,  and  Dm  that 
for  the  model.  While  they  do  not  give  expressions  similar  (5.9),  Brown  and 
Gangi  in  their  published  article  (1963)  do  note  that  scaling  keeps  the 
number  of  wavelengths  between  scaled  distances  constant,  which  is  equivalent 
to  equation  (5.6.c),  i.e. , loss-less  wavelengths. 


In  practice,  there  are  severe  restrictions  on  using  these  relations  in 
a simulator  model.  Paraphrasing  Iizuka  (1968) : 

(a)  the  values  of  er  and  loss  tangent  for  the  model  are  limited 

(b)  the  dimension  scale  factor  has  to  be  such  that  the  model  is  not 
too  small  or  too  bulky 

(c)  ranges  of  er  and  loss-tangent  by  mixing  materials  is  limited 

(d)  changing  loss-tangent  or  er  is  complex,  it  being  possible  to 
minimize  one  while  changing  the  other 

(e)  materials  must  not  be  noxious,  explosive  or  too  costly. 


5 . 3 laboratory  lithospheric  models  and  similitude 

5.3.1  Idealized  waveguide  with  electrical  constants  uniform,  with  depth 

Brown  and  Gangi  (1962,  1963)  and  Gangi  (1966)  made  laboratory  measure- 
ments with  a simulated  waveguide,  the  simulation  representing  an  "average" 
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set  of  conductivity  values  deduced  from  their  proposed  profiles  (See 
Chapter  2) . 

The  model  of  the  lithosphere  to  be  simulated  is  shown  in  the  Figure  5.1. 

In  their  measurements,  the  earth  model  was  assumed  to  have  a depth  De  = 35  km, 
with  er2  = £re  = 4.  Since  their  conductivity  depth  model  (Chapter  2)  showed 
a conductivity  <?2  ~ CTm  that  might  vary  with  frequency,  appropriate  scaling 
ratios  were  employed. 

An  aluminum  (or  carbon)  plate  represented  the  overburden  or  upper  crust 

and  an  aluminum  plate  was  used  for  the  lower  crustal  boundary  (Moho) . The 

crustal  waveguide  was  a salt  water  solution  whose  conductivity  could  be 

5 

varied.  A frequency  scale  factor  of  fm/fe  = 10  was  chosen  for  simplicity 

in  the  models.  The  value  of  em  for  water  was  80.  Since  waveguide  distances 

of  1000  km  were  desired  to  be  simulated,  and  a conveniently-sized  tank  was 

desired  for  the  laboratory,  a distance  (or  depth  Dm/De)  scale  factor  works 
5 

out  to  be  4.47  x 10  ; thus  if  De  were  35  km,  a tank  depth  of  7.82  cm  or 
3.08  inches  results.  Distances  were  plotted  in  units  of  Hg  or 
(»e  " De'  % = Dm). 

5. 3. 1.1  Sample  results 

In  Figure  5.2  is  shown  measurements  on  the  model  at  300  MHz,  in  the 

tank  3.08  inches  in  depth,  of  the  vertical  (relative)  electric  field  strength 

over  a range  of  2H  to  17H  (when  H is  or  He,  the  conductivity  cm  being 

0.07  Si/m.  The  loss  tangent  pm  = ,0525  which  is  small;  the  corresponding 

/.  -6 

plane  wave  attenuation  is  then  = 6Oir0mA  /e^  =*  3.30  x 10  n/m  = 0.28 
dB/km,  These  values  of  loss-tangent  for  the  earth  pe  and  ag  work  out  to  be 
the  same,  since  they  are  scaled. 


4 
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In  comparison  with  theory,  Figure  5.2  shows  bracketing  calculations 

as  dots  or  crosses;  the  dots  correspond  to  0m  = .05  Si/m  for  which  the 

loss  tangent  p_  is  0.0375  for  f = 300  MHz  just  as  it  is  for  the  scaled 
in  m 

-8 

value  of  ag  = 2.5  x 10  Si/m  at  a frequency  of  3KHz.  The  crosses  are 

for  calculations  for  a = 0.10  Si/m  for  which  the  loss-tangent  p = 0.075 

m ^ ^m 

—8 

at  300  MHz,  corresponding  to  the  earth  values  = 5 x 10  Si/m  at  3 KHz. 

The  calculations  bracket  well  the  measured  curve. 

Brown  and  Gangi  (1962,  1963)  and  Gangi  (1966)  show  the  effects  on 
field  distributions  when  simple  inhomogeneities  were  introduced.  These 
included  (a)  a step  barrier  in  the  lower  boundary  to  correspond  to  the 
shallower  depth  of  the  Moho  under  the  sea  at  the  continent-sea  boundary; 

(b)  a bar  in  the  upper  boundary  to  represent  a geologic  fault  zone,  and 

(c)  a cylinder  to  correspond  to  a lens  or  intrusion  of  mountains  into  the 
granite.  While  idealized,  the  results  showed  potential  interference  pat- 
terns compared  with  those  when  no  discontinuity  was  present.  (See  the 
original  papers  for  detailed  results) . It  should  be  noted  that  the  ocean 
crust  may  differ  markedly  in  electrical  properties  from  that  of  the  continen- 
tal crust  (Cox,  1971,  1978). 

5.3.2  Lithospheric  waveguide  with  constants  varying  with  depth  (inhomogeneous) 

An  inhomogeneous  lithospheric  waveguide,  i.e. , a waveguide  whose  elec- 
trical parameters  vary  with  depth  has  been  investigated  theoretically  in 
the  past  seven  or  eight  years  but  relatively  little  has  been  reported  on 
simulation  experiments.  While  the  depth  variations  may  be  represented  by 
thin  horizontal  slabs,  it  is  a bit  more  satisfying  to  simulate  the  smooth 
(sic)  depth  profiles  by  suitable  materials  whose  properties  may  be  smoothly 


changed.  About  the  only  extended  model  investigations  reported  are  those 
of  Iizuka  (1967,  1968,  1969,  1971,  1972).  I have  found  no  theoretical 
comparisons  with  his  measurements.  His  results  extend  earlier  ones  on 
antenna  impedance  with  stratified  films  of  Iizuka  and  King  (1962  a,b) . 

The  simulator  chamber  used  an  agar-agar  solution  (gel)  into  which  a 
solution  of  sodium  chloride  (NaCl)  could  diffuse  to  change  the  concentration. 
Two  chambers  were  built,  one  for  antenna  impedance  measurements  and  a second 
larger  one  for  field  strength  measurements.  The  technique  of  fabricating 
inhomogeneous  media  was  described  in  Iizuka* s 1967  paper;  a summary  of  the 
work  on  impedance  and  field  strength  appears  in  his  1972  paper. 

The  variation  of  loss  tangent  and  relative  dielectric  constant  as  a 
function  of  frequency  are  shown  in  Figure  5.3  for  various  media  such  as 
agar-agar,  agar-agar  with  ethyl  alcohol,  tap  water,  0.1  mol  NaCl  solution 
and  ethyl  alcohol.  For  agar-agar  NaCl  affects  loss  tangent;  ethyl  alcohol 
affects  er. 

Agar-agar  has  a high  dielectric  constant  and  for  simulation  use  it 
can  reduce  the  dimensions  of  the  simulation  chamber.  For  various  reasons 
(See  Iizuka  1968,  Chapter  III),  a modelling  frequency  of  114  MHz  was  chosen. 
The  properties  of  agar-agar  with  concentration  of  NaCl  at  114  MHz  are  shown 
in  Figure  5.4.  On  the  Figure  equations  for  least  squares  fit  are  given  for 
er  and  tan  6,  as  a function  of  the  concentration  C of  the  NaCl.  The  agar- 
agar  solution  was  1.3%  agar-agar  and  98.7  tap  water.  The  concentration  of 
NaCl  was  the  weight  in  grams  in  1000  cu  cm  of  agar-agar.  If  this  concen- 
tration C is  less  than  3,  then  the  approximations  simplify  to: 
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er  - 78.30  + 0.073  C (5.14) 

p = tan  5 =>  0.1569  + 0.4633  C (5.15) 

for  114  MHz. 

A choice  of  fm  can  be  made  by  consideration  of  Figure  5.3,  and  equation 
(5.6)  or  (5.13).  From  Figure  5.3.  er  of  agar-agar  is  near  78  for  frequencies 
less  than  4000  MHz,  then  decreases  at  higher  frequencies.  The  loss-tangent 
starts  to  increase  near  4000  MHz.  This  limits  the  agar-agar  choice  of  n fm, 
n being  the  dimension  scale  factor  re/rm,  and  the  original  full  scale  fre- 
quency fe.  Thus  from  equation  (5.13) : 


The  relative  dielectric  constants  of  the  agar-agar  em  vs  fm  are  known  from 
Fiqure  5.3.  In  Figure  5.5,  the  model  simulation  frequency  fm  is  plotted 
vertically  against  the  product  of  the  full  scale  frequency  fe  and  the  dimen- 
sional scale  factor  n for  various  values  of  relative  dielectric  constant  Ere 
of  the  full  scale  system.  The  lowest  loss-tangent  of  agar-agar  is  that  get 
with  no  NaCl.  It  is  plotted  on  the  left  of  Figure  5.5  vs  the  model  frequency 
^m*  These  curves  can  determine  model  parameters.  If  the  lithosphere  wave- 
guide is  to  be  operated  at  fe=5KHz,  with  ere  * 4,  and  if  the  distance  factor 
5 

is  n * 10  , then  n fe  = 500  (if  fe  is  in  Mhz) . The  model  frequency  then  is 
fm  = 114  MHz  at  which  the  minimum  loss-tangent  is  0.13. 

Various  profiles  of  the  lithosphere  itself  were  deduced  from  values  pub- 
lished up  to  about  1965  and  combined  into  one  plot  from  Iizuka  (1968)  and 
are  shown  in  Chapter  2.  A typical  simulated  profile  is  shown  in  Figure  5.6 
for  a frequency  fe  of  5 KHz  for  loss  tangent,  the  relative  dielectric  constant 
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c''f)  simulation  of  typical  lithospheric  profile  with  thot  in  motel  a ar-cvar  chamber  (From  Iizuka.1969  ) 
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ere  being  4 and  constant  with  depth.  The  depth  of  the  lithosphere  Dm  was 
30  km,  and  a modelling  scaling  factor  of  10^  was  used  so  that  the  simulation 
chamber  was  30  cm  deep,  Dm  = 0.3  meters.  An  aluminum  plane  represented  the 
large  loss-tangent  upper  crust,  and  a water  solution  of  0.2  mol  NaCl  repre- 
sented the  mantle.  The  simulated  profile  of  the  conducting  layer  was 
obtained  by  diffusion  of  the  NaCl,  through  a porous'  film  into  the  agar-agar 
at  the  bottom. 

Iizuka  (1968,  Chapter  XX)  describes  how  to  fabricate  profiles  of 
various  shapes  including  (a)  complimentary  error  function,  (b)  error  function, 

(c)  linear  profile,  (d)  periodic  profile,  and  (e)  miscellaneous  profiles. 

5. 3. 2.1  Field  strength  measurements 

It  should  be  mentioned  that  the  experimental  modelling  chamber  of  Iizuka 
could  be  arranged  also  with  non-parallel  top  and  bottom,  a vertical  step 
profile  (similar  to  Brown  and  Gangi,  1963) , a vertical  slant  profile,  or 
with  simultaneous  horizontal  and  vertical  profiles. 

5. 3. 2. 1.1  Non-parallel  plate  region  with  homogeneous  medium 

The  geometry  is  shown  in  Figure  5.7.  Table  I shows  the  dimensions  and 
electrical  properties  of  the  medium  for  various  geometers.  Note  that  for 
loss  tangents  of  0.036,  0.23  and  0.88,  the  corresponding  conductivities  are 
0.073,  0.47  and  1.78  Si/m,  the  latter  approaching  that  for  sea  water. 

The  resulting  electric  field  strength  E in  Figure  5.8  is  shown  plotted 
as  E R vs  the  range  R for  a quarter-wave  monopole  radiator.  Note  maxima 
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and  minima  for  the  lowest  loss-tangent.  The  smaller  the  angle  the  smaller 
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Dimensions  and  electrical  properties  of  the  medium  in  a non-parallel  plate 
region. 
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FIG.5B  PLOT  OF  E^R  vs.  R IN  UNPARALLEL  PLATE  REGION 
FILLED  WITH  A CONDUCTING  MEDIUM  AT  114  MHz. 
SEE' THE  PARAMETERS  IN  TABLE  5.1  (1.2 IWJ 
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the  attenuation.  As  conductivity  increases,  the  attenuation  becomes  less 
sensitive  to  angle  changes.  That  is,  at  high  conductivity  the  propagation 
depends  more  on  the  losses  of  the  medium  than  by  the  geometry  of  the  bottom. 
This  is  useful  to  know  for  operation  at  sloping  shore  lines. 

5. 3. 2. 1.2  Field  strength  measurements  in  an  inhomogeneous  medium 

For  the  lithospheric  guide  model  an  agar-agar  gel  represented  the 
lithosphere,  covered  with  a metal  plate  top.  A quarter-wave  monopole 
extended  through  the  plane  into  the  medium.  As  such,  the  configuration 
was  similar  to  a parallel  waveguide  with  metallic  walls  except  that  the 
bottom  was  replaced  by  an  inhomogeneously  conducting  wall.  Figure  5.9 
is  a representative  result.  The  depth  profile  of  er  and  loss-tangent  is 
shown  in  Figure  5.9  (b) , while  the  amplitude  and  phase  of  the  magnetic 
field  strength  are  shown  in  Figure  5.9  (a).  Note  the  model  was  0.383 
X thick,  yet  the  field  strength  decreases  monotonically  with  depth.  The 
field  was  measured  at  R = 50  cm;  a perfect  guide,  less  than  half  a wave- 
length thick  would  have  a field  magnitude  invariant  with  depth. 

5. 3. 2. 1.3  Comment  - Images 

Figure  5.9  represents  measurements  on  a model  lithosphere  for  which 
the  upper  wall  was  a metal  plate  above  the  inhomogeneous  lithospheric  model. 
In  essence,  this  could  be  considered  one-half  of  a lithosphere  having  a 
symmetrical  inhomogeneous  profile,  the  upper  half  being  the  image  of  the 
lower  one  shown.  Since  the  excitation  was  by  a vertical  monopo’e,  the 
result  would  be  that  of  a center-driven  dipole,  giving  maximum  field  strength 
at  the  center  where  the  conductivity  (loss  tangent)  is  minimum. 
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Several  of  the  more  recently  postulated  profiles,  such  as  those  of 
Keller  or  of  Housely  discussed  in  Chapter  2,  can  be  approximated  by  quasi- 
symmetrical  profiles.  Calculations  of  Field  and  Dore  (1973)  show  maximum 
field  strengths  near  that  for  minimum  conductivity. 

Advantage  could  be  taken  of  this  image  effect  along  with  the  various 
inhomogeneous  model  profiles  achieved  by  Iizuka  using  the  agar-agar  gel 
with  sodium  chloride  and  ethyl  alcohol  diffusion  discussed  in  Section  5.3.2. 
The  uniform  profiles  of  Brown  and  Gangi  (1963)  with  saline  water  could 
likewise  be  treated  - the  scaling  parameters  would  have  to  be  reconsidered 
See  Chapter  6. 

The  use  of  images  for  measuring  wall  losses  is  also  discussed  in 
Chapter  6.  Here,  however,  the  question  of  similitude  scaling  application 
is  difficult  to  answer  because  the  walls  are  necessarily  water  instead  of 
metal , in  order  to  make  distance  or  depth  dependent  measurements  in  the 
"walls"  be  the  fields  excited  by  dipoles  in  the  lithosphere  or  in  the 
walls.  Nevertheless,  the  principles  will  be  demonstrated  and  the  simpler 
measurements  understood  before  proceeding  step  by  step  to  more  complex 
cases. 

This  leads  to  a discussion  of  modelling  materials. 

5.4  Modelling  Materials 

5.4.1  Typical  High-conductivity  Materials 

Fritschkneckt  (1971)  gives  a table  of  typical  high  conductivity 
materials  which  might  be  considered  for  modelling  purposes,  and  is  shown 
in  Table  5.2.  Listed  are  values  for  metals,  other  solids,  and  some 
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Material 

Conductivity 
mho  m 

Metals 

Aluminum 

1.44  ' 107 

Aluminum,  commercial  alloys 

l.22-3.o0-l0 

Brass 

1.2- 1.3-1  O' 

Bronze,  phosphor 

0.97- 1 .04-  107 

Bronze.  “Everdur* 

3.4-5.3-10* 

Copper 

3.41  - 10’ 

Copper,  tube  (DHP  copper* 

5.01-  107 

Lead 

4.55  ■ 10* 

Magnesium 

2.17  - 10' 

M.ingantn 

2.27  - 10’ 

Mercury 

1.04  ■ 10* 

'■'tool,  nonmagnetic  stumlcss 

1.35-2.14 ' io' 

Zinc 

1 .02  • 10* 

Other  solids 

Carbon 

2.0- 2.4  • 104 

Graphite 

0.7-  1.2  - 105 

Kpow.  commercial  silver  loaded 

5.0-  10* 

Epovy.  commercial  copper  loaded 

2.5  • 104 

-Lav-up*  resin,  graphite  loaded 

0-  1.0  • '.0* 

’•Vast,  bronze  loaded 

-1.25-  \04 

Silicone  rubber,  silver  loaded 

1.3  - 10* 

L t qu  ids 

\mmomum  chloride. 

l'*%  by  weight  in  water  it  14*C 

17.4 

tapper  sulfate,  is  above 

J.  J 

llvdmohlonc  acid,  as  above 

02.0 

•vadium  chlorine,  ss  above 

12.0 

Sulfuric  acid,  as  above 

30.2 

Sea  a. iter 

1-  *> 

YAA/f  Of  SOLoCNT 


.BENZENE 


-CRMULA  '0f£j^CamjCrviTY\S0LtjaL' 


y\oxKlA'LMmmaLlT''  1^^'! 


2 296  532'0 


*R0  ; FLAM 


K:HnWI.I.H 


CUTE 

53/GAL  , 

no  a 

0XYLSN6  iCj* 


p-CYMENE CH3C«H4CHiCH 


NAPHTHALENE  'C 


! TTPAHYCRONAPHTHAL£Ne!c1o  h. 


PfNENE  iC.0H. 


OIL  OF  TURPENTINE 


flam 

3i3/  CO  9 

334  0 j 

Bgaam 


p-CRESOL 

CHjC*ruOH 

5. 

0IETHYL  ether 


01  ;SO  amyl  £TH£R 


' ICH^CMlCHjIjuO  2 91 


139 


'anisoyl  chloride 

CH,OC,h4. 

F 

ACUTE  | 

i 

170  39 

j-CRESYL  . 


m-CRESYL  METHYL 


o-CRESYL  METHYL 


OlOXANE 


METHYLAL 


n-auTYL  ACETATE 


,0-0- c. 


iCH}C,h«)zO  3 57 


ICH.C, 


(CHjC«H,),0  4 03 


OCH 


Ch,;0CH,1 


9 I 


2 • ;0" 


CHjCOOICHjIjCHj  5 CO  ! 2 55  • '0 


flam  I 62/2509  '182  0 


'98  9 : 


FLAM.EXPL  I 14/100  9 1 101  3 


CHRO  I FLAM 


acute  : 1540/1004  '118  IT 


l CHRO 


12  93/2309  '26  3 


30  -AMYL  aoQP-ONATE  CH,CH;C00;CHg,CK|  4 23 


'SO-AMYL  n-PuTYRATE  ;CH.,,;CH,J.COO  3 95 


SO-amyl  SO-yALERiANATE  'CH,)*iCH,],CHjCOO|  3 52  | 


ME’hyl  BENZOATE 


5 03 


8ENZYL  3ENZ0ATE  OCOCCH.O  ,4  83 


S6/1OO9 

113  5 1 

unstable 

CHLOROFORM  Oh  0 7) 


CARBON*E~RAChlORiCE  CCI. 


3FNTACHL0RE  HANE  ChC»,  CCI, 


C:S  acetylene  OIChlCRiCSCm  CJCh,;/ 


M 

, ACUTE 

100/1  ,b 

61  26  : 

CYCLOHEXANE 


mE"hylcY 


OE  : AHYHRCNAPH  ^UA_EN£ 


• H . HARD 
J ■ GCCC 

F , FAIR 


*S«SOLUBlL.TY 

C'CCST 


'A3L2  =.3  low  iialactric  constant  solvents  (iron  lizuka, 1 -6l ) 


i 




- 5.26  - 

t 

liquids.  When  large  scaling  ratios  are  required  for  models,  a difficulty 
appears.  One  would  like  to  have  a material  with  conductivity  between  that 
for  saline  water  and  carbon,  for  certain  frequency  ranges.  For  distance- 
dependent  field  measurements,  liquids  or  gels  (such  as  agar-agar  gels  of 
Iizuka)  are  indicated. 

5.4.2  Low  Dielectric  Constant  Solvents 

While  searching  for  materials  of  low  dielectric  constants,  Iizuka 
(1961)  tabulated  the  results  of  his  survey,  shown  here  as  Table  5.3. 

Dioxane  was  indicated  but  was  too  toxic  and  dangerous,  so  other  solvents 
were  used.  The  materials  were  desired  in  obtaining  the  performance  of 
antennas  in  stratified  media.  The  solvents  were  separated  by  thin  films 
such  as  nylon  (having  low  dielectric  constants  and  conductivities) . The 
horizontal  stratification  effects  are  reported  by  Iizuka  (1961)  and  Iizuka 
and  King  (1962b). 

5.4.3  Higher  Dielectric  Constant  Liquids 

For  scaled  modelling  purposes, . higher  dielectric  constant  liquids  are 

desired.  For  this  purpose,  water  is  most  often  used,  or  ethyl  alcohol  or 

a mixture  of  the  two,  providing  a range  of  20  to  80  or  so  for  relative 

dielectric  constant.  By  adding  salt,  the  conductivity  can  then  be  varied 
-4 

from  10  to  about  10  Si/m. 

While  less  desirable  for  showing  smooth  profiles  than  the  agar-agar 
mixtures  of  Iizuka,  inhomogeneous  profiles  cam  be  simulated  by  horizontally 
stratified  layers  of  solvents  with  differing  electrical  constants  according 
to  the  profile  desired.  Such  stratification  might  ultimately  be  considered 
in  the  schemes  of  Chapter  6 as  one  of  the  steps  in  building  up  from  simpler 


models 
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6.  Theory  and  Experiment  In  lithospheric  Propagation:  Models 
6.1  Introduction 

The  possibility  of  wave-guide-like  propagation  of  electromagnetic  waves 
in  the  earth's  lithosphere  has  been  studied  extensively  by  numerous  investi- 
gators. The  theoretical  models  are  all  quite  similar  and  are  based  on  ideali- 
zations of  cross  sections  of  the  earth's  crust  such  as  that  given  by  Levin 
[1966]  and  shown  in  Fig.  6.1.  These  idealizations  in  general  include  (a)  a 
planar  slab  of  a dielectric  material  with  low  conductivity  to  represent  a 
lithosphere  composed  of  granitic  rock  and  basalt,  (b)  an  upper  layer  of  quite 
highly  conducting  material  corresponding  to  the  overburden  on  the  continent 
and  the  ocean  basin  beyond  this,  and  (c)  a lower  also  highly  conducting  layer 
to  represent  the  mantle.  Refinements  sometimes  added  include  the  air  and  the  ion- 
osphere above  the  overburden  and  ocean  as  shown  in  Fig.  6.2,  but  their  presence 
usually  has  a negligible  effect.  The  excitation  of  electromagnetic  waves  in 
the  dielectric  slab  and  their  reception  at  distant  points  has  taken  several 
forms.  The  first  and  most  commonly  described  source  consists  of  a vertical 
antenna  in  a borehole  extending  into  the  lithosphere  or,  short  of  this,  ending 
in  the  overburden  at  some  distance  from  the  lithosphere.  Wait's  analysis  [1966] 
shows  that  if  the  vertical  antenna  is  not  in  the  lithosphere  or  very  close  to 
it,  adequate  excitation  is  not  possible.  He  suggests  [1966,  p.923]  that  "ever/ 
effort  should  be  made  to  locate  both  the  (vertical)  source  dipole  and  the  obser- 
ver within  the  waveguide.  Otherwise,  some  other  mechanism  of  excitation  should 
be  used  such  as  a horizontal  electric  dipole."  This  suggestion  is  followed  by 
Frieman  and  Kroll  [1973]  whose  Investigations  include  excitation  not  only  by  a 
vertical  electric  dipole  but  also  by  a horizontal  electric  dipole  at  the  upper 
boundary  of  the  lithosphere.  They  conclude  that  "the  horizontal  electric  dipole 
case  proves  to  be  more  interesting",  that  it  "permits  the  utilization  of  much 
larger  antennas  and  hence  more  efficient  antenna  design."  They  suggest  that  "a 
horizontal  electric  dipole  is  far  superior  to  a vertical  electric  dipole." 


o . - 


Fig. 6.1  Profile  of  earth's  crust  according  to  S.B,  Levin 
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The  method  of  analysis  of  Frieman  and  Kroll  is  similar  to  that  of  Wait 
and  others  in  that  the  relevant  waveguide  modes  are  determined.  Their  study  of 
the  propagating  characteristics  of  the  several  modes  led  them  to  the  conclusion 
that  for  a horizontal  antenna  and  receiver  on  the  ocean  floor  the  behavior  of 
the  field  "should  be  very  similar  to  that  of  horizontal-dipole  excitation  near 
the  plane  boundary  of  two  semi -inf inite  conducting  media."  However,  it  is  well 
known  [BaHos,  1966;  King  and  Sandler,  1977]  that  transmission  along  such  an 
interface  is  primarily  by  so-called  lateral  waves  for  which  a complete  theory  is 
available.  To  date  it  has  not  been  applied  to  an  interface  between  two  half- 
spaces with  the  properties  of  sea  water  and  basalt  or  of  the  earth’s  overburden 
and  granitic  rock.  It  appears,  therefore,  that  for  the  idealized  theoretical 
model  consisting  of  a dielectric  slab  between  conducting  regions,  analytical 
and  numerical  solutions  of  Maxwell's  equations  are  available  and  should  suffice 
to  determine  the  nature  of  propagation  between  a pair  of  vertical  or  horizontal 
dipoles  in  the  slab  or  near  its  boundary  in  the  upper  conducting  layer.  In 
particular,  it  should  be  possible  to  calculate  the  electric  field  at  the  re- 
ceiver as  a function  of  distance,  frequency,  and  the  assumed  electrical  parameters 
of  the  media.  Alternatives  to  such  computations  include, on  the  one  hand,  the 
very  difficult  and  costly  full-scale  experiments  and,  on  the  other  hand,  suitably 
scaled  and  constructed  experiments  that  satisfy  all  of  the  conditions  of  electro- 
magnetic similitude.  Model  measurements  have  the  advantage  over  theoretical 
computations  that  they  permit  a variety  of  realistic  modifications  of  the  over- 
idealized theoretical  models. 

In  the  absence  of  magnetic  materials  the  conditions  for  modelling  are: 
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where  f is  the  frequency,  r is  a characteristic  length,  e - e0er  is  the  PerTIlit' 
tivity  and  a is  the  conductivity.  The  subscript  m denotes  the  model,  the  sub- 
script e the  earth.  Unfortunately  such  modelling  is  severely  limited  by  the 

distances  involved  and  by  the  materials  and  frequencies  available  for  the 
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simulation  and  the  experiments.  A model  of  practical  size  requires  a scale 

factor  in  the  distances  of  the  order  of  r /r  ■ 10^  . With  e - 4 to  8 for 

e m er 

the  lithosphere  and  e ■ 80  for  the  water  solutions  most  commonly  used  in  the 

1/2 

model,  20  to  10.  It  follows  that  - (£mr/£er)  (rm/re}  " 

4 x 10~5  to  3 x 10-^  . Since  accurate  electrical  measurements  are  difficult 

unless  f = 1 GHz,  it  follows  that  f « (45  to  30)  KHz.  That  is,  completely 
m e 

scaled  modeling  of  the  lithosphere  and  its  upper  and  lower  bounding  regions  is 

limited  to  actual  frequencies  that  satisfy  f » 45  KHz.  While  this  condition 

presents  no  practical  limitation  on  VLF  propagation  in  the  lithosphere  using 

vertical  dipole  excitation,  it  excludes  all  of  the  high  frequencies  which  may 

be  the  most  promising  for  lateral  wave  transmission  using  horisontal  dipole 

excitation.  The  second  condition  for  similitude  requires  that  for  eer  ■ 4,  a&/am  m 

2 2 e. 

fr/fr  »2x  10  If  the  lithosphere  is  modeled  with  water  solutions  as 
m m e e * 

-4 

was  done  by  Brown  and  Gangi  and  by  Iizuka,  * 10  to  10  Si/m,  so  that 

o - 2 x 10"1Q  to  2 x 10“5  Si/m.  This  includes  the  entire  useful  range.  On  the 
e 

-2  -3 

other  hand  for  an  overburden  with  a * 10  to  10  Si/m,  a m 3000  to  300  Si/m 

e m 

for  which  convenient  materials  are  difficult  to  find.  The  same  is  true  for 

the  ocean  with  a - 4 Si/m  and  o = 4 x 10^  Si/m,  with  e__  * 80* 

e m er 

Actual  models  constructed  by  Brown  and  Gangi  [1963]  were  designed  to  simu- 
late a frequency  f * 2 to  4 KHz.  The  lithosphere  was  quite  accurately  scaled 

by  salt  water,  the  overburden  or  ocean  and  the  mantle  less  accurately  by  metal 
or  carbon  plates.  For  waveguide  propagation  when  the  transmitting  and  receiving 
antennas  are  both  vertical  and  embedded  in  the  lithosphere,  a less  than  exact 
scaling  of  the  properties  of  the  upper  and  lower  conducting  media  may  be  adequate 
so  long  as  the  principal  attenuation  is  due  to  losses  in  the  lithosphere  and  not 
in  the  bounding  walls.  Note,  however,  that  the  use  of  solid  plates  for  the  over- 
burden or  ocean  is  neither  convenient  nor  sufficiently  accurate  when  the  exci- 
tation is  by  a vertical  or  horizontal  dipole  located  at  some  point  above  the 
surface  of  the  lithosphere. 
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The  interesting  models  developed  by  Iizuka  [1971,1972]  have  similar 
general  limiations  as  those  of  Brown  and  Gangi  with  respect  to  the  frequency 
range  and  the  properties  of  materials.  They  are  also  restricted  to  excitation 
by  dipoles  in  the  model  lithosphere.  However,  they  do  provide  the  means  for 
producing  a specified  non-uniform  conductivity  profile  across  the  lithosphere 
and  for  providing  boundaries  that  are  not  plane  parallel.  In  a restricted 
range  of  frequencies  and  with  excitation  by  a vertical  antenna  within  the 
lithosphere  Iizuka' s models  offer  a powerful  tool  for  obtaining  information 
about  wave  propagation  in  the  earth  insofar  as  the  assumed  profiles  are  good 
approximations  of  reality.  Unfortunaely  they  are  not  readily  generalized 
to  permit  the  use  of  vertical  or  horizontal  dipoles  in  the  model  overburden  or 
ocean  and  for  determining  the  electromagnetic  field  in  this. 

6.2  Suggested  Experimental  Studies  in  Lithospheric  Propagation 

Measurements  on  models  combined  with  the  principle  of  similitude  provide 
a powerful  and  practical  alternative  to  full-sized  experiments  in  the  earth 
and  to  theoretical  studies.  Owing  to  the  very  large  scale  factors  involved 
and  the  limitations  of  experimental  techniques  and  available  materials  the 
method  is  for  all  practical  purposes  limited  to  low  frequency  applications 
in  the  actual  earth.  Even  for  these  the  accurate  modelling  of  the  lithosphere 
virtually  precludes  a corresponding  accuracy  in  the  modelling  of  the  overburden 
or  the  ocean  and  so  requires  the  exciting  and  receiving  antennas  to  be  in  the 
model  lithosphere.  In  the  actual  earth  this  is  possible  only  in  boreholes 
through  the  overburden  or  into  the  ocean  bed  with  vertical  dipoles.  In  view  of 
the  theoretically  derived  conclusion  by  Frieman  and  Kroll  [1973]  that  horizontal 
dipoles  are  a very  promising  alternative,  future  experimental  studies  should 
be  concerned  with  a comparative  investigation  of  the  characteristics  of  both 
horizontal  and  vertical  electric  dipoles  when  located  in  the  ocean  near  its 
interface  with  the  lithosphere  and  not  only  with  vertical  dipoles  in  the  litho- 
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sphere.  Suggested  experiments  might  include:a)  vertical  transmitting  and 
receiving  dipoles  in  the  model  lithosphere  and  at  increasing  heights  in  the 
model  overburden  or  ocean;  b)  a vertical  transmitting  antenna  in  the  model 
lithosphere  and  a horizontal  receiving  dipole  in  the  model  overburden  or  ocean 
near  the  interface  with  the  lithosphere;  and  c)  horizontal  transmitting  and 
and  receiving  antennas  in  the  model  overburden  or  ocean.  In  each  case  the 
receiving  dipole  is  to  be  moved  to  measure  the  magnitude  and  phase  of  the 
electric  field  as  the  distance  from  the  transmitting  antenna  is  increased. 

Further  measurements  are  indicated  to  determine  the  effect  of  non-parallel 
boundaries  of  the  model  lithosphere,  boundary  surfaces  between  the  lithosphere 
and  the  overburden  or  ocean  that  are  not  smooth  but  undulate  or  are  otherwise 
irregular  or  discontinuous,  a layer  of  properly  modelled  sediment  on  the 
ocean  bottom  including  the  effect  of  embedding  both  transmitting  and  receiving 
antennas  at  various  depths. 

Since  communication  between  two  horizontal  dipoles  in  the  overburden  or 
ocean  near  its  interface  with  the  lithosphere  is  primarily  by  lateral  waves 
[Banos,  1966;  King  and  Sandler,  1977],  experiments  that  involve  such  antennas 
must  be  designed  to  take  maximum  advantage  of  the  somewhat  unusual  character- 
istics of  this  type  of  transmission.  In  particular,  note  must  be  taken  of  the 
fact  that  when  the  horizontal  transmitting  antenna  is  not  far  from  the  inter- 
face, the  most  useful  component  of  the  electric  field  is  the  radial  one  in  a 
direction  along  the  axis  of  the  dipole.  Furthermore,  the  optimum  frequency 
for  transmission  is  well  above  the  low  frequency  range.  Since  accurate  model- 
ling according  to  the  principle  of  similitude  cannot  be  carried  out  above  full  scale 
frequencies  of  two  to  three  kilohertz,  it  is  not  practical  to  study  lateral  wave 
transmission  in  a comprehensive  manner  by  models  that  satisfy  the  equations  of 
similitude  as  formulated  for  the  lithosphere.  While  this  means  that  a completely 
scaled  set  of  experiments  is  unavailable,  it  does  not  mean  that  meaningful  and 
illuminating  information  cannot  be  obtained  from  suitably  designed  model  experi- 


ments.  These  must,  of  course,  be  planned  and  interpretted  in  conjunction  with 
the  theory  of  lateral-wave  transmission. 

In  carrying  out  the  suggested  experiments  it  may  be  convenient  and  expedi- 
ent to  separate  the  transmission  losses  in  the  model  lithosphere  from  those  in 
the  model  overburden,  ocean,  and  mantle  . This  can  be  done  initially  by  making  use 
of  a perfect  dielectric  (e.g.  air  or  styrofoam)  in  the  model  lithosphere  and  sea 
water  for  the  bounding  regions.  Subsequently  a dielectric  with  losses  and 
specified  conductivity  profiles  comparable  to  those  believed  to  obtain  in  the 
earth  can  be  introduced  in  the  simpler  and  more  crucial  experiments. 

6.3  Some  Possible  Laboratory  Models  for  Lithospheric  Experiments 

The  simple  model  used  by  Brown  and  Gangi  [1963]  and  the  much  more  complicated 
one  developed  by  Iizuka  [1971,1972]  stress  the  accurate  modelling  of  the  litho- 
sphere but  represent  the  overburden  or  ocean  and  the  mantle  by  highly  conducting 
metal  or  carbon  plates.  In  order  to  study  the  significance  of  the  bounding 
regions  in  their  effect  on  the  transmission  and  simultaneously  to  make  possible 
the  embedding  and  moving  of  antennas  in  them,  it  is  possible  to  model  them  with 
a liquid  like  water  or  alcohol  or  a mixture  of  these  to  obtain  relative  permit- 
tivities between  20  and  80  and  to  add  salt  to  vary  the  conductivity  in  the 
-4 

range  from  10  to  10  Si/m.  The  model  lithosphere  can  be  composed  initially  of 
air  or  styrofoam  with  negligible  losses  and  subsequently  by  wood  or  alcohol  with 
higher  conductivity.  The  bounding  styrofoam  can  be  smooth,  undulating,  or 
arbitrarily  irregular.  Above  the  air  or  styrofoam  is  a metal  plate  that  can  be 
parallel  to  the  surface  of  the  water  or  inclined  at  an  arbitrary  angle.  It 
serves  as  an  image  plane  for  the  entire  structure  below  it.  This  effectively 
provides  a model  lithosphere  that  has  twice  the  thickness  between  the  metal 
plate  and  the  surface  of  the  water.  It  is  bounded  above  and  below  by  salt 
water  to  model  both  the  overburden  or  ocean  and  the  mantle. 

A sketch  of  the  suggested  model  is  shown  in  Fig.  6.3.  The  vertical 
transmitting  and  receiving  monopoles  when  combined  with  their  images  act  as 


-6.8- 


transmitting 

monopole 


receiving 

metal  plate  monopole  i 

'I 
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(a)  Schematic  diagram  of  apparatus 
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(salt  water) 
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model  overburden  or  ocean 
(salt  water) 


(b)  Schematic  diagram  of  apparatus  with  image 


Fig.  6.3  Suggested  apparatus  for  lithospheric  and  lateral-wave 
transmission. 
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dipoles  in  the  center  of  the  model  lithosphere.  The  receiving  monopole  travels 
in  a slot  to  measure  the  amplitude  and  phase  of  the  vertical  electric  field 
from  near  the  transmitter  to  distant  points.  The  upper  plate  can  be  inclined 
and  the  water-styrofoam  surface  can  be  made  undulating  ot  otherwise  irregular, Fig. 6. 

Also  shown  in  Fig.  6.3  are  transmitting  and  receiving  dipoles  in  the  salt 
water.  In  the  figure  they  are  horizontal  but  in  practice  they  can  be  oriented 
to  be  vertical  or  at  any  desired  angle.  The  receiving  element  can  be  used  to 
measure  the  amplitude  and  phase  of  the  electric  field  in  the  model  overburden 
or  ocean.  All  components  of  the  electric  field  can  be  determined.  Since 
the  dipoles  have  identical  images  which  are  symmetrically  located  in  the  model 
mantle,  the  field  actually  measured  is  due  to  both  the  source  and  the  image. 

However,  the  contributions  are  readily  separated  theoretically,  so  that  the 
field  due  to  one  source  alone  can  be  studied. 

In  order  to  carry  out  the  suggested  measurements  in  a tank  approximately 
3x2  meters,  a frequency  in  the  range  from  0.3  to  1 GHz  is  required.  This 
gives  between  three  and  ten  wavelengths  in  air  for  studying  the  decay  of  the 
electric  field  in  the  model  lithosphere  when  excited  by  a horizontal  dipole 
in  the  overburden  or  ocean  or  by  a vertical  monopole  in  the  model  lithosphere. 

With  the  former  the  effect  of  varying  the  distance  between  the  antenna  and  the 
bounding  surface  can  also  be  investigated.  The  general  arrangement  shown  in 
Fig.  6.3  provides  the  means  for  carrying  out  the  suggested  experiments  and  so 
obtaining  direct  information  on  the  transmission  of  electromagnetic  waves 
through  the  lithosphere  and  of  lateral  waves  along  its  boundary  under  a vari- 
ety of  circumstances  arranged  to  simulate  possible  conditions  in  the  earth. 
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Such  a program  of  experimentation  and  its  correlation  with  theory  would 
seem  to  be  a logical  and  necessarv  continuation  and  extension  of  the  basic 
work  by  Brora  and  Gangi  and  especially  by  Iizuka. 
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7.  Summary  and  Conclusions 

Propagation  in  the  earth's  crustal  waveguide  has  been  studied  by  numerous 
investigators.  The  properties  of  available  conceptual  depth  profiles  of  elec- 
tric constants,  mainly  conductivity,  are  reviewed.  The  propagating  region  was 
assumed  to  be  a planar  slab,  bounded  on  top  by  an  ocean  or  overburden  and  on 
the  bottom  by  the  mantle  or  Moho,  the  bounding  regions  being  of  very  large 
loss- tangent. 

The  waveguide  theories  extant  are  then  summarized.  Earlier  analyses 
assumed  a slab  of  uniform  depth  profile.  Many  of  the  investigators  illustrate 
the  modal  properties  for  modes  of  low  attenuation,  others  give  samples  of 
field  calculations.  Usually  the  exciting  sources  were  vertical  electric 
dipoles  (VED)  extending  through  bore-holes  in  the  overburden  (or  ocean)  into 
the  lower  loss-tangent  crustal  waveguide.  Analyses  for  non-homogeneous  pro- 
files included  the  modal  solutions  of  Wait  for  an  exponential  conductivity 
profile,  and  those  of  Schwering,  Peterson,  and  Levin  for  an  inverse-square 
depth  profile  of  both  dielectric  constant  and  conductivity.  Full  wave  modal 
solutions  were  employed  by  Field  and  Dore  using  recent  conceptual  ”bell-shaped" 
profiles  of  Keller  and  of  Housley  which  showed  minima  and  very  much  smaller 
conductivities  at  depths  of  the  order  of  10  km. 

Peking  VED  excitation  in  the  waveguide,  Wait  suggested  HED  excitation, 
a suggestion  which  was  followed  up  by  Frieman  and  Kroli  who  showed  that  for 
ocean  emplaced  antennas  at  the  ocean-lithosphere  interface,  the  HED  was, 
superior  to  the  VED.  Excitation  should  be  like  that  for  HED  excitation  near 
the  plane  boundary  of  two  semi-infinite  media.  For  dipoles  in  the  ocean 


near  the  ocean-air  boundary,  King  and  Sandler  extended  the  work  of  Banos. 


Propagation  is  by  means  of  lateral  waves.  This  work  was  extended  to  the  case 
where  the  air  was  replaced  by  the  lithosphere.  Calculations  were  made  for  a 
wide  range  of  parameters  for  the  radial  component  of  electric  field  from  the 
HED,  which  component  is  the  useful  one  especially  for  intermediate  and  large 
ranges.  The  mechanism  may  be  described  as  "down-under-and-up"  in  consonance 
with  the  popular  notation  of  "up-over-and-down"  for  HED  sources  in  the  sea 
near  the  sea-air  boundary.  The  advantages  and  disadvantages  of  the  lateral 
wave  transmission  were  discussed  for  the  ocean-lithosphere  case. 

The  relations  for  the  scaling  of  models  according  to  electrodynamic 
similitude  are  presented.  Examples  are  given  from  the  work  of  Brown  and  Gangi 
for  homogeneous  depth  profiles  where  the  modelled  crustal  slab  (saline)  water 
is  bounded  on  top  by  metal  plates.  For  non-homogeneous  profiles,  examples 
are  given  from  the  work  of  Iizuka  who  used  an  aaar-agar  gel  with  diffusion  of 
sodium  chloride  to  vary  the  loss-tangent  with  depth.  The  top  boundary  was  a 
metal  plate,  the  bottom  being  the  higher  conductivity  sodium  chloride  diffuser. 
Measurements  were  made  with  planar  and  sloping  boundaries. 

To  study  the  effects  of  bounding  regions  on  transmission  including  losses 
and  simultaneously  make  measurements  of  the  effects  of  embedding  and  moving 
antennas  in  them,  the  laboratory  models  should  have  "walls"  composed  of  a 
liquid  like  water  or  alcohol  or  a mixture  of  them  which  can  vary  the  relative 
dielectric  constant  from  20  to  80  and  to  add  salt  to  vary  the  conductivity 
from  10  to  10  Si/m.  Above  such  a wall  would  be  the  lithosphere,  initially 
air  or  styrofoam,  and  subsequently  a substance  like  wood  or  alcohol  to  in- 
crease the  complex  dielectric  constant.  Above  the  "lithosphere"  is  a metal 
plate  which  can  be  parallel  or  inclined.  The  styrofoam  boundary  can  be 


undulating  or  ir regular.  The  metal  plate,  which  has  movable  probes,  acts 
as  an  image  plate. 

It  is  recommended  that  an  orderly  build-up  of  modelling  experiments  be 
initiated  with  a laboratory  model  just  described,  and  continuing  with  more 
complicated  depth  profiles.  Although  available  materials  may  limit  accurate 
modelling  according  to  the  principles  of  similitude,  meaningful  and  illumi- 
nating information  can  be  obtained  for  a better  understanding  of  lithospheric 
transmission  including  lateral  waves  and  losses. 

The  more  complicated  crustal  depth  profiles  require  a much  more  elaborate 
set-up  possibly  patterned  after  that  of  Iituka.  For  measurements  and  their 
interpretation  then,  there  will  be  a need  for  scientists  of  considerable 
experience  including  support  personnel.  The  program  should  be  long  range 
in  anticipation  of  the  development  of  more  accurate  depth  profiles  from  full 
scale  geophysical  experiments  and  analyses. 
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